Mott physics and first-order transition between two metals in the normal state phase 

diagram of the two-dimensional Hubbard model 
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For doped two-dimensional Mott insulators in their normal state, the challenge is to understand 
the evolution from a conventional metal at high doping to a strongly correlated metal near the Mott 
insulator at zero doping. To this end, we solve the cellular dynamical mean-field equations for the 
two-dimensional Hubbard model using a plaquette as the reference quantum impurity model and 
continuous-time quantum Monte Carlo method as impurity solver. The normal-state phase diagram 
as a function of interaction strength U , temperature T, and filling n shows that, upon increasing 
n towards the Mott insulator, there is a surface of first-order transition between two metals at 
nonzero doping. That surface ends at a finite temperature critical line originating at the half-filled 
Mott critical point. Associated with this transition, there is a maximum in scattering rate as well 
as thermodynamic signatures. These findings suggest a new scenario for the normal-state phase 
diagram of the high temperature superconductors. The criticality surmised in these systems can 
originate not from a T=0 quantum critical point, nor from the proximity of a long-range ordered 
phase, but from a low temperature transition between two types of metals at finite doping. The 
influence of Mott physics therefore extends well beyond half-filling. 

PACS numbers: 71.27.-|-a,71.30.-|-h 
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I. INTRODUCTION 

The discovery of quantum oscillations in layered high- 
temperature superconductorai"— has brought renewed at- 
tention to the normal phase below the superconducting 
dome. The challenge posed is to understand how these 
systems evolve as a function of carrier concentration in 
the absence of the superconducting phased. At large 
doping, the copper-oxide layers are metallic. On the 
other hand, at zero doping these systems are Mott in- 
sulatorsSiS, i.e. systems that are insulating due to the 
strong screened Coulomb interactions between the elec- 
trons. Thus, as a function of the carrier concentration, a 
Mott metal-insulator transition (MIT) occurs, i.e. a tran- 
sition driven by the on-site Coulomb repulsion between 
the electronsiSiii. We are challenged to understand the 
Mott transition driven by carrier concentration in the 
normal phasei^. 

At the theoretical level, to address this physics we need 
to study the competition between the kinetic band ef- 
fects, that delocalize the electrons in the lattice, and the 
effects due to screened interactions that localize them. 
These two ingredients are present in the Hubbard model, 
which contains a kinetic energy term that describes the 
band structure, plus an interaction term that represents 
the screened local Coulomb repulsion between electrons 
occupying the same site of the lattice. More than fifty 
years of intense research shows that this simple model has 
a surprisingly rich set of solutions that may help capture 
the complex behavior observed in nature with minimal 
assumptions. 



Even the simple-looking Hubbard model poses a 
formidable challenge to theory. Dynamical mean-field 
theory (DMFT) has played a central role in advancing 
our knowledge of the physics contained in the Hubbard 
model and has proven to be a useful method to inves- 
tigate the Mott metal insulator transition. In essence, 
dynamical-mean-field theory embeds a single sit o^'^d^ 
or a cluster^'5"^^ in an infinite self-consistent bath of 
non-interacting electrons. The latter problem is then 
amenable to accurate numerical solutions. This approach 
can also be justified from a variational perspective^. 
DMFT with a single site immersed in a bath is exact 
in infinite dimensio n^'^d"^ and provides a quite accurate 
mean field description of materials^^ and of the Mott 
transitioniS in three dimensions. The picture is that the 
Mott insulator can be driven towards a correlated metal- 
lic state through a first-order transition by tuning corre- 
lation strength, or temperature, or dopingi^i^Sr— . Here 
the fundamental unifying concept is the first-order transi- 
tion. That transition dominates the entire normal phase 
diagram of the model, and is relevant at finite tempera- 
ture even if the actual Mott transition is hidden by the 
onset of some long-range order. 

Something similar occurs in the Fermi liquid descrip- 
tion of the normal state of weakly correlated materials. 
There, the zero-temperature Fermi liquid fixed point can 
be masked by a long-range ordered state that is the true 
ground state. Nevertheless, as long as one is not too 
close to the transition to long-range order, the normal 
state is well described by the finite temperature Fermi liq- 
uid that emerges from the zero-temperature Fermi liquid 
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fixed point. While the Fermi liquid provides a mean-field 
description of a phase where the self-energy is analytical 
(7mS oc w^), DMFT also allows a phase, such as the 
Mott insulator, where the self-energy is not an analytical 
function of frequency^!. 

For strongly-interacting layered materials, such as the 
above-mentioned cuprate high temperature supercon- 
ductors, or for the layered organic conductors, the de- 
scription provided by the single-site DMFT solution of 
strongly correlated electron models is valid only at very 
large doping or at temperatures that are in fact much 
too large to be interesting. This is because single-site 
DMFT neglects spatial fluctuations, and for these sys- 
tems the short-range correlations play an important role 
due to their low dimensionality. For example, the ob- 
served pressure-driven first-order Mott transition in or- 
ganic superconductors of the BEDT family22r— has the 
wrong slope in the pressure-temperature plane when 
compared with the single-site DMFT solution. This is 
corrected^iH by cluster extensions of DMFTi^^iii^i-^ 
that incorporate short-range magnetic correlations in the 
theoretical description. Similarly, quantum Monte Carlo 
calculations on the square lattice at half-filling show that 
a Mott gap opens up at finite temperature independently 
of antiferromagnetism starting around^ U ~ 6t, with U 
the interaction strength and t nearest-neighbor hopping. 
It is only with extensions of DMFT on a plaquette that 
this result is recovered^"—. We stress that as larger and 
larger cluster sizes are considered^, the description of 
the normal state will become more and more accurate at 
lower temperature, but it will also fail when long-range 
order sets in. Hence, we argue that calculations, such as 
ours, that are based on a plaquette and reproduce the 
U ~ 6t result mentioned above are an accurate mean- 
field description of the normal sate. Hence, this is the 
type of approach that we adopt in the present paper. 

Cluster extensions of DMFT are however computation- 
ally more expensive than their single-site counterpart. 
Nevertheless, useful physical insights on the phase dia- 
gram of the Hubbard model have been obtained within 
that framework and have helped elucidate some im- 
portant properties of the copper oxide superconductors. 
Among the main results, we mention pseudogap forma- 
tion at low doping^"—, the presence of antiferromag- 
netism and superconductivity^!^"—, and the momen- 
tum dependence of the electronic properties along the 
Fermi surface^i^"— . We caution that cluster extensions 
of DMFT come in two varieties: The Dynamical Clus- 
ter Approximatioii^i (DCA) and cellular DMFT^i^. 
For simplicity we did not distinguish between both ap- 
proaches in the previous discussions. They generally give 
qualitatively similar results. Here we shall employ cellu- 
lar DMFT. 

Returning to the normal state, the overall picture for 
the Mott transition in cluster extensions of DMFT have 
been lacking because cluster methods, as already men- 
tioned, are computationally expensive. However, a few 
years ago there was a breakthrough. The continuous- 



time quantum Monte Carlo technique in the hybridiza- 
tion expansion^"— now provides us with the possibility 
to explore the full phase diagram of the model within the 
cluster DMFT method. This is because this algorithm 
considerably reduces the computational time, allowing 
both unprecedented accuracy and studies over a wider 
range of parameters. Recent investigations using this 
method and other continuous-time quantum Monte Carlo 
techniques^, have determined the interaction driven 
Mott MIT, revealing sharp modifications to the single- 
site picture^i^i^. The transition driven by carrier con- 
centration, more relevant for the high-temperature su- 
perconductors, is also currently under intense investi- 
gation^MLi^i^iH-^. Motivated by the physics of the 
cuprates, almost all the studies focus on the large in- 
teraction regime where the Mott gap is well developed 
and they consider the effects of different band structure 
parameters in order to capture the striking particle-hole 
asymmetry observed in those compounds. This direction 
of research produced several important results, however 
a unified scenario for the Mott transition realized in clus- 
ter extensions of DMFT for the Hubbard model was still 
missing. 

In our recent work"^- and in the present study, we take 
a different approach to address this problem. Instead of 
focusing on a region of the parameter space where the 
interaction strength is large, we map out the whole nor- 
mal state phase diagram of the two-dimensional Hubbard 
model as a function of interaction, temperature and dop- 
ing. To capture the general features of the localization- 
delocalization Mott physics, we simply consider nearest- 
neighbor hopping. 

Our strategy is motivated by one of the key lessons of 
the single-site DMFT solution of the Hubbard model: a 
first-order transition governs the entire phase diagram of 
the model. We ask whether the Mott transition still dom- 
inates the phase diagram and how it is modified by the 
short-range correlations taken into account through cel- 
lular DMFT. As mentioned above, in the Hubbard model 
the Mott transition can be tuned by temperature, by the 
interaction strength (interaction driven MIT, relevant for 
the layered organics) and by the carrier concentration 
(doping driven MIT, relevant for the cuprates). Fun- 
damentally, one is exploring the same strong coupling 
physics from different perspectives. At half-filling, the 
cluster DMFT solutions of the model do display a first- 
order transition, as a function of the interaction strength, 
between a correlated metallic state and a Mott insula- 
tor^-ia. 

What is the fate of this first-order transition as a func- 
tion of the carrier concentration? As reported in our 
recent work*"" and in great detail here, the full mapping 
of the normal-state phase diagram reveals that upon in- 
creasing the carrier concentration towards the Mott in- 
sulator, there is indeed a surface of first-order transition. 
This is distinct from the surface of first-order transition 
found at fixed [/ as a function of second-neighbor hop- 
ping in earlier work^®'^^. Thus, as in the single-site case. 
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in cellular DMFT the finite temperature Mott transition 
provides the key to understand the whole phase diagram. 
In addition, contrary both to single-site DMFT and to 
conventional wisdom, the first-order transition, emerg- 
ing at the half- filled Mott critical endpoint, progressively 
moves to large doping as the interaction strength in- 
creases. Hence the transition can occur between two 
metallic states, one of which evolves continuously from 
the Mott insulator, a surprising effect brought about by 
the short-range correlations. 

The first-order transition and the associated critical 
line that we find are not only new, they also provide a 
unified picture for a host of previously known results. 
For example, the first-order transition was not detected 
at large values of the interaction by previous works most 
likely because the critical line moves rapidly to lower 
temperature with increasing interaction strength, falling 
below the temperatures accessible to date. Hence, pre- 
viously reported features of the phase diagram of the 
Hubbard mode l^^'^^i^*^" — , whether it is thermodynamic 
properties, scattering rate, momentum differentiation or 
other, appear in a different light. As we shall see, our 
analysis allows one to identify these features as precur- 
sors of the first-order transition, hence calling for their 
re-evaluation. It is in this decisive sense that our work 
goes beyond -and is compatible with- previous studies. 

Compared to our recent letter™, the present work ex- 
tends and provides further results on the thermodynam- 
ics, on the spin susceptibility, on results for other values 
of the parameters and many additional details and a re- 
fined interpretation. In Section |TT] we briefly discuss the 
model and method. The peculiar normal-state phase dia- 
gram of the two-dimensional Hubbard model is the topic 
of Section IIIII In Section IIVI we discuss the first-order 
transition that controls the physics of the phase diagram. 
Section |V] contains the thermodynamic properties of the 
phases separated by the transition. The characterization 
of these phases continues in Section I VII where a signa- 
ture of critical behavior through a large scattering rate 
is found near the transition. We identify the physical 
origin of this critical behavior in the Section IVIIl In Sec- 
tion |VllT] we reexamine the phase diagram and propose a 
possible new scenario for the surmised criticality in high- 
temperature superconductors. It can originate not from 
a quantum critical point, nor from the proximity of a 
long-range ordered phase, but from a finite-temperature 
transition at finite doping coming from the influence of 
Mott physics well beyond half-filling. 



II. MODEL AND METHOD 

The Hamiltonian for the two-dimensional Hubbard 
model on a square lattice reads 



Here Ci^ and cf^ operators annihilate and create elec- 
trons on site i with spin a, and Uia- — cf^Cia is the 
number operator. In this article we focus on the physics 
arising from the proximity to a Mott insulating state, 
so we use a simple hopping amplitude tij between near- 
est neighbors only, so that the bare dispersion, obtained 
from the Fourier transform of tij, is given by e(fc) = 
—2t[cos{kx) + cos{ky)]. U is the energy cost of double 
occupation at each site of the lattice, is the chemical 
potential which, for a given value of U, controls the oc- 
cupation n — ^/N "^^^{rii^) where N is the number of 
sites. 

We solve this model using cellular dynamical mean- 
field theoryi^ii^i^. This approach maps the lattice prob- 
lem Eq. ([I} onto a quantum impurity problem consisting 
of a cluster of sites, here a 2 x 2 plaquette, embedded 
in a bath that is determined self-consistently in such a 
way that infinite lattice and plaquette have the same self- 
energy. The action of the quantum impurity problem 
(plaquette coupled to the bath) is given hy^ 
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(1) 



s = s. 



dr dr'^Vl?(T)A«fl,(T,T')^i?'(T'), 

Jo Jo j^j^, 

(2) 

where Sc is the action of the cluster, R labels the cluster 
sites, and A is the 4x4 bath hybridization matrix. The 
self-consistency condition that fixes the bath hybridiza- 
tion matrix A is 



A(ia;„) =ia;„ + fj, - ic - Sc(«w„) 

1 



E 



iujn + t{k) - Sc(jw„) 



(3) 



where Sc is the cluster self-energy matrix, tc is the pla- 
quette hopping matrix, i(fc) is the lattice matrix of hop- 
ping in the supercell notation and fc runs over the reduced 
Brillouin zone of the superlattice. 

There are a variety of techniques to solve the quan- 
tum impurity problem Eq. . Here we use the re- 
cently developed continuous-time quantum Monte Carlo 
method^"—. This approach relies on the Monte Carlo 
summation of all diagrams generated by the expansion of 
the cluster-bath (impurity) action Eq. ^ with respect to 
the hybridization A. Ref. [6ll contains the details of the 
present implementation. 

The CTQMC method is a powerful and perfectly 
adapted technique to map out the phase diagram of the 
Hubbard model, because it permits to readily access all 
regions of doping, from large doping down to the Mott 
insulating state, the intermediate to low temperature 
regime that have so far proven inaccessible with usual 
Hirsch-Fye quantum Monte Carlc*"^ for instance, and a 
large range of \J with high efficiency^. In addition, this 
method is statistically exact and does not have errors 
.associated with time discretization or bath parametriza- 
tion, and therefore we can determine the phase diagram 
of the Hubbard model with unparalleled accuracy. 
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To obtain high quahty data, required to determine the 
phases of the system, we typically perform 5 x 10^ Monte 
Carlo sweeps per processor, averaged over 40 processors. 
When necessary, for example close to phase boundaries, 
we may do up to 10^ sweeps per processor averaged on 
64 processors. Cellular DMFT is an iterative approach, 
and in generic regions of parameter space we have stud- 
ied the cluster observables converge in less then 20 — 30 
iterations, but hundreds may be necessary close to phase 
boundaries. 

In the actual CTQMC implementation, it is useful to 
write cluster quantities, such as the hybridization A that 
appears in Eq. ([2|), the cluster self-energy, or the clus- 
ter Green function, in a diagonal form. Since cellular 
DMFT applies open boundary conditions on the cluster, 
the cluster momentum K is not in general a good quan- 
tum number. However, in a 2 x 2 plaquette, one still has 
GA^ symmetry and the one-dimensional irreducible repre- 
sentations of the subgroup C2^, equivalently noted by K, 
can be used. Then the cluster quantities take a diagonal 
form with indices (0,0), corresponding to the represen- 
tation Ai, (tt, tt) corresponding to the representation A2 
and (tt, 0), (0, tt) corresponding to the representations Bi 
or B2. For example, the bath hybridization function in 
cluster momentum basis reads: 
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^(7r,7r 



(4) 



The diagonal elements are the eigenvalues of the matrix 
containing the on site, nearest neighbor, and next-nearest 
neighbor cluster quantities in real space. Physically they 
can be thought as an average of the lattice quantity over 
a coarse grained approximation of the Brillouin zone. In 
analogy with multiband models, we refer to these clus- 
ter momenta K as "orbitals". When exploiting this in- 
terpretation, one should keep in mind, however, that in 
the cellular DMFT scheme the cluster momenta K are 
coupled through both the on-site interaction U and the 
hopping matrix t, which in cluster momentum basis takes 
the form: 



i{k) - 



tl2 il3 

^22 ti3 

^13 ^33 ti2 

^13 ^12 ^44 



(5) 



where tn = —t (2 + cos{2kx) + cos{2ky)), ^22 ~ 
t{cos{2kx) — cos{2ky)), ^33 — —t{cos{2kx) — cos{2kyj), 
^44 = t{2 + cos(2fca;) -I- cos(2fcy)), ti2 = itam{2kx), 
ti3 = itsin{2ky), and kx,ky S [— 7r/2, 7r/2]. For refer- 
ence, the non-interacting local density of states on the 
orbitals obtained from 



G{uj)-- 
is shown in Fig. [1] 



dk a; + /.t — t{k) 



(6) 



— K=(0,0) 

- - K=(jt,0) 




FIG. 1: Non interacting local density of states pk(!-^) ~ 
-l/-KlmGK{i^) of the orbitals K = (0,0), (7r,0), (tt, tt) (solid, 
dashed and dot-dashed lines respectively). 



CTQMC is an imaginary time technique, thus for the 
interacting case the local density of states is not di- 
rectly accessible and has to be computed by analytical 
continuation schemes that rely on additional approxima- 
tions. Nevertheless, the characterization of the phase 
diagram of the normal phase of the Hubbard model - 
which is the topic of this study - can unequivocally be 
determined from observables directly measured in the 
CTQMC method, like the particle density n, the cluster 
Green's function and self-energy in Matsubara frequency. 
Here we restrict our study to that kind of quantities. 
From a methodological viewpoint, this is a necessary step 
before using other methods to extract additional informa- 
tion. 



III. PHASE DIAGRAM 

In this section we describe the normal-state phase dia- 
gram of the two-dimensional Hubbard model determined 
in detail by cellular DMFT calculations^^. 

The parameter space is three dimensional and consists 
of all possible values of the interaction strength [/, the 
temperature T and the chemical potential /i (or, equiv- 
alently, the doping 5 defined as 5 = 1 — n). To iden- 
tify the different phases and locate their boundaries or 
crossover lines, one must scan the phase diagram with 
an adapted fine grid. Despite the algorithmic break- 
throughs described in the previous section, that have be- 
come available only in the last few years, it took several 
million of CPU hours to be able to determine the whole 
normal-state phase diagram of the two-dimensional Hub- 
bard model. 

Specifically, we carried out scans at constant values of 
U across the Mott transition, varying the chemical po- 
tential for several temperatures. Because the system is 
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FIG. 2: Chemical potential /i, interaction U, temperature T phase diagram of the two-dimensional Hubbard model obtained by 
cellular DMFT. Because of particle-hole symmetry it is symmetric with respect to the fj, = plane. Cross-sections at constant 
U are shown. Dark-gray (blue) shaded regions represent the coexistence of two phases. Light-gray (yellow) areas denote the 
onset of the Mott insulator state (MI), characterized by a plateau in the occupation at n = 1. When these two regions overlap, 
a metal- insulator transition takes place (different shade of gray). Otherwise, the coexistence regions occur between two different 
metals. Projections on T = 0, and ^ = planes are also shown (full lines and dashed lines respectively). Open dots mark the 
extrapolated T = values of Ud and Uc2- A critical line Ta (dotted line) originates at the half-filled Mott critical endpoint 
^''mit (full dot) and moves to progressively low temperatures and high doping as U increases. 



particle-hole symmetric {f is set to zero) , we limited our 
study to hole doping only {fi < 0, S < 1). For each 
value of U, we performed calculations in a wide doping 
range, between and roughly 25%. We found that scan- 
ning doping in intervals sometimes smaller as 0.002 was 
necessary to identify the sequence of phases that occur 
upon doping the Mott insulator. As far as temperature 
is concerned, we mostly focused on the low temperature 
regime in the decade 1/100 < T/t < 1/10. 

We summarize our results in the temperature T versus 
interaction strength U and chemical potential fj, phase 
diagram shown in Fig. [2] The scans at constant U that 
we use to explore the three dimensional parameter space 
are visible as cross sections in the phase diagram and 
correspond to the distinct T — /i planes presented iuFigH 

The figures show two colored (shaded) regions in cross 
sections at constant U in ([/, T, /i) space. The first re- 
gion (light gray/yellow) corresponds to the onset of the 



Mott insulating phase, characterized by a plateau in the 
occupation at n = 1. Just outside this region there is a 
metallic state. The second region (dark gray /blue) is the 
portion of parameter space where two different phases 
coexist. By interpolation between the distinct planes at 
constant U, one can therefore obtain the volume of both 
the Mott phase and the phase coexistence region. 

It is important to understand that what we mean by 
coexistence region is a region where, in addition to the 
thermodynamically stable phase, a metastable phase can 
exist. Since we do not have access to the free energy, 
we cannot determine precisely the first order boundary, 
so instead we draw the region where metastable states 
exist. In practice, if we change chemical potential slowly 
from small to large values or vice- versa, we can end up 
in different phases, one of which is thermodynamically 
unstable. 

Going back to our phase diagram, let us first concen- 
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(e) (f) (g) (h) 



FIG. 3: Phase diagram in the temperature T versus chemical potential ^ plane for diflterent values of the interaction strength 
U. Each panel corresponds to a cross-section at constant U of the previous figure. As before, light gray (yellow) area represents 
the Mott insulating phase (MI). Dark grey (blue) area represents the coexistence of two phases and is bounded by the spinodal 
lines jj,ci{T) and ^c2{T). A star symbol marks the end of the coexistence region at finite temperature; by following this point 
as a function of U gives rise to the dotted line Tcr of Fig. O Other symbols denote some of the points actually computed in our 
study. Three phases can be distinguished: overdoped phase (open circles), Mott insulator (crosses), underdoped phase (filled 
squares). The former is metallic. In the latter some regions of the Brillouin zone are gapped and others gapless, and the system 
is compressible. 



trate on the Mott phase. Two basic features emerge. 
First, as expected on physical grounds, a threshold value 
for the correlation strength, called J7ci, is a prerequisite 
for the emergence of Mott insulating state. We estimate 
this point at {U,T,fi) w (5.3^,0,0). Second, the projec- 
tions on T = plane (marked as light gray /yellow line on 
the 3D phase diagram) reveal the characteristic V shape 
onset of the Mott insulating state, where the tip of the V 
shape is Ud (only half of the V shape is shown in Fig. [2] 
because of particle-hole symmetry). 

Our most surprising results come, however, from the 
full mapping of the region (dark gray /blue) where two 
phases coexist in the {U,fi,T) space. We begin with 
a description of the peculiar shape of this region. Let 
first consider the fi = plane, where the model is half 
filled. In that plane, the coexistence region has a triangu- 
lar shape that bends as U is increased. The boundaries 
are the spinodal lines Uci{T) (thick dashed blue line), 
where the metastable insulating phase disappears, and 
Uc2{T) (thin dashed blue line), where the metastable 
metallic phase ceases to exists. These spinodals end 
at a critical value of C/, called J/mit, with coordinates 
(t/, T, /i) « (5.95, 0.08, 0). In the ^ = plane, the region 
of coexistence marks a portion of parameter space where 
the Mott insulator coexists with a metal. 

Surprisingly, this is in general not so in the 3D phase 
diagram where the coexistence region is not fully in- 
cluded in the region where the Mott phase exists. To 
show this, let us consider scans at constant U . The re- 



gion where two phases coexist naturally extends from the 
/i = plane to finite values of the chemical potential /i 
and is enclosed by the spinodal surfaces iici{U,T) and 
fic2{U,T) that end at a critical line Tcr (dotted blue line 
in Fig. [5]). That line is the continuation of the critical 
point iU,T,fj,) w (5.95,0.08,0) from half filling to finite 
doping. The dark grey (blue) shaded areas in the phase 
diagrams correspond to the intersection between the vol- 
ume delimited by these surfaces and the cross sections at 
constant U. 

Let us look in more detail at the coexistence regions. 
In Fig. [3] we show constant U cuts of Fig. [51 For values of 
the interaction strength U larger than C/mit, the region of 
coexisting phases displays again a triangular shape bend- 
ing towards the Mott insulator (i.e. towards low doping). 
As we discuss in the next section, two metallic phases co- 
exist in that region. The coexistence begins at the tem- 
perature Tcr, marked by a star symbol in Fig. [3l and is 
delimited by the spinodals /Xci(T) and fic2{T). This is 
apparent from the scans at U/t = 6.0, 6.2. For larger U, 
the coexistence region narrows and drops to temperatures 
lower than our current numerical capabilities. On the 
other hand, as the interaction strength U decreases to- 
wards Uun, the triangular region grows and approaches 
the 11 — plane and its summit at the top reaches that 
plane at C/mit- As U is further decreased in the range 
Uc2{T = 0) < U < C/mit, the coexistence regions show a 
rainbow-like shape (cf scan at C/ = 5.8t) with the /ici(T) 
line on the outer part of the arc and the fic2 (T) line on 
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the inner part of the arc. Further lowering U in the re- 
gion Uci{T = 0) < U < Uc2{T = 0), only the spinodal 
/ici exists until Uci{T — 0) where it vanishes. 

Our results can be summarized as follows. First, within 
the volume in (U — fi—T) parameter space enclosed by the 
spinodal surfaces ^ici{T, U) and Hc2(T, U), one can obtain 
two solutions of the cellular DMFT equations. Therefore 
a first-order transition surface has to occur in the blue 
region. The precise location of the first order surface can 
be computed by comparing the free energy of the two 
solutions, a task beyond the present work. 

Second, the two solutions merge at a finite temperature 
Tcr where the spinodal surfaces fici{T, U) and Hc2{T, U) 
end. This implies that a second-order transition line 
(dotted blue line in Fig. ^ exists at the end of the sur- 
face of first-order transition. Except for the critical Mott 
endpoint J7mit, the second-order critical line takes place 
outside the Mott insulating phase, so it occurs for finite 
values of doping. 

Third, our most striking finding is the nature of the 
phases that are found to coexist. The first-order tran- 
sition surface originates at the half-filled Mott critical 
endpoint and moves progressively away from half-filling 
as the interaction U increases. Therefore the first-order 
transition can occur between a metallic phase and a Mott 
insulating phase, or, unexpectedly, between two metallic 
phases, one of which evolves continuously from the Mott 
insulator. The former case (metal to insulator transi- 
tion) occurs whenever the volume of the Mott insulating 
phase intercepts the volume of the phase coexistence in 
the ([/, /x,T) parameter space. As discussed above, this 
occurs at the /i = plane (i.e. at half filling), where 
the metallic and the insulating solutions of the cellular 
DMFT equations coexist. In the T — fj, planes too, the 
Mott insulator can coexist with a metal, as indicated 
in Figs. I2I3I by different shade of gray (corresponding 
to the light gray/dark gray - yellow/blue overlap) for 
U/t = 5.6, 5.8. Even though in some region of the phase 
diagram the first-order transition occurs between a metal 
and a Mott insulator, this is generally not the case. For 
sufficiently large U (U > C/mit), there is a first-order 
transition between two metals, as can be seen in Fig. [2] 
where for U > Umit there is no intercept between the co- 
existence region (dark gray /blue) and Mott region (light 
gray /yellow) . In all cases shown in the phase diagram, 
however, the fj,c2iU,T) spinodal surface denotes the van- 
ishing of a metallic solution, while the fici{U^ T) spinodal 
surface does not necessarily coincide with the end of Mott 
insulating phase. 

In the following, we refer to the metallic phases sep- 
arated by the first-order transition as underdoped (UD) 
and overdoped (OD) phase. Here these terms indicate 
whether the level of doping of the phase is below or above 
the level of doping at which the first-order transition oc- 
curs. The underdoped phase is closest to the Mott insu- 
lator. The connection of these names to the maximum of 
the superconducting dome will be discussed in Sec. IVIIII 



IV. FIRST-ORDER TRANSITION AT FINITE 
DOPING 

The most important results in the normal-state phase 
diagram surveyed so far is the first-order transition sur- 
face in the 3D parameter space (U, T, fj,) with the different 
phases it separates. The focus of this section is a system- 
atic analysis of this first-order transition. We begin with 
the numerical results that demonstrate the existence of 
such a transition, and then turn to the thermodynamic 
behavior of several observables near its boundary. 

Fig. m demonstrates the existence of a finite T first- 
order transition. We plot the occupation n as a function 
of the chemical potential for a wide range of tempera- 
ture and eight values oi U. It is important to stress that 
each curve n{p) is obtained by keeping all model param- 
eters fixed except /i. A plateau in the curves at n(/i) = 1 
appears above the critical coupling Uci{T) and signals 
the incompressible Mott insulating phase. The marked 
temperature dependence of the n{fj,) curves in the UD 
phase, as opposed to the OD phase, demonstrates that 
the effects of strong correlations manifest themselves at 
low energy in the UD phase. 

For U/t — 5.6,5.8,6.0,6.2 and low temperatures, the 
occupation clearly shows hysteretic behavior as a func- 
tion of fi, an unequivocal fingerprint of the first-order 
nature of the transition. To find the hysteresis loop, 
we use an iterative procedure on a fine grid of fi points, 
where a converged solution is used as seed for the next 
value of /i. We obtain the lower branch of the hysteresis 
cycle by increasing /i starting from large enough dop- 
ing. This branch corresponds to the metallic branch 
of the overdoped phase and is indicated by full sym- 
bols in Fig. 01 For U/t ~ 5.6, which lies in the region 
Uci{T = 0) < [/ < t/c2(T = 0), this branch continues 
up to half filling, at = 0. In contrast, for U larger 
than Uc2{T = 0), the metallic branch of the OD phase 
shows an upward jump at finite doping. The OD metal- 
lic branch endpoint defines fj,c2- On the other hand, to 
compute the upper branch of the hysteresis loop, we de- 
crease /i starting from the Mott insulating solution at 
half filling. This branch is indicated by open symbols 
in Fig. m As /X varies in the Mott plateau, this branch 
remains constant at n = 1, then it evolves continuously 
into a compressible metal (i.e. the underdoped phase), 
and subsequently undergoes a downward jump at finite 
doping. The endpoint of the upper branch defines /id. 

The hysteresis region delimits a region of space pa- 
rameters {U — T — n) where two solutions of the cellular 
DMFT equations coexist. The crossing of the free energy 
of the two solutions in this volume bounded by /Zci {T, U) 
and /ic2 {T, U) defines the first-order transition surface. 
Determining this surface and finding out which is the 
most stable phase in the coexistence region necessitates 
to compute the free energy. This task requires a large 
investment in computing time and is beyond the scope 
of the present work. 

The continuous evolution from the Mott plateau to the 
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FIG. 4: Occupation n versus ^ for several values of interaction strength U . The data shown are for temperatures T /t = 1/10 
(blue triangles), 1/25 (green squares), 1/50 (red circles), 1/100 (black diamonds). When two solutions are found to coexist, the 
solutions obtained following the metallic and the insulating solution are indicated as full and open symbols respectively, pci 
marks the vanishing of the underdoped phase. /ic2 signals the disappearance of the overdoped phase. The arrows indicate the 
values of ^ci and ^c2 at the lowest temperature shown in the panel. The plateau in the occupation at n = 1 signals the onset 
of the incompressible Mott state. Note that ^ci in general occurs at finite values of doping. The inset in the panel for U = 6.2t 
shows the n(/i) curves for temperatures T/t — 1/40 (orange down triangles), 1/50 (red circles), 1/60 (magenta left triangles), 
1/64 (violet right triangles), which are above the second order critical temperature Tcr ~ 1/65* at which the two phases merge. 
Note the sigmoidal shape of the n{n) curves as Tcr is approached from above. 



compressible UD phase seems to occur towards a discon- 
tinuous change in the compressibility (first derivative of 
n{fi)) at T = which would correspond to a second order 
transition, or quantum critical point. This is apparent es- 
pecially at large U. At finite temperature, only higher 
order derivatives are perhaps discontinuous. 

The fj,c2{T,U) surface always coincide with the van- 
ishing of the metallic state (namely, the OD phase), as 
revealed by the jump of the lower branch at finite dop- 
ing. In sharp contrast, the ptci{T, U) surface does not in 
general mark the vanishing of the Mott insulating state: 
as shown in Fig. |4l the sudden jump at /Xd occurs at 
finite doping! These two observations lead to the con- 
clusion that the first-order transition can actually occur 
between two types of metallic phases. Since these phases 
have the same symmetry, the first-order surface can end 
at a critical line at finite temperature Tcr (dotted blue line 
in the 3D phase diagram of Fig. [5]). As Tcr is approached 
from below, the hysteresis loop narrows and vanishes at 
Tcr- On the other hand, as Tcr is approached from above, 
the curves show a sigmoidal profile, a clear signature of 
proximity to a second-order critical point. The inset of 



Fig. 4(e) reveals this remarkable feature. 



Hysteresis in the occupation curves n{^)^ found for 
U/t = 5.6,5.8,6.0,6.2 and low temperatures, provides 
a direct demonstration of the first-order nature of the 
transition. The curves at U/t = 7.0, 9.0, 12.0 do not show 
hysteresis down to T/t = 1/100, but several signatures 
indicate that the system is getting close to the transition 
and hence suggest that the hysteresis loop falls below the 
lowest temperature we investigate. To reach this conclu- 
sion, first we notice that a close inspection of the n{fi) 
curves at U/t = 6.0,6.2, which lie above the critical in- 
teraction Umit ~ 5.95t, shows that the temperature Tcr 
dramatically decreases as U increases. Indeed, the crit- 
ical line begins at the Mott point at T/t ^ 0.08 and 
decreases to T/t ^ 0.015 in going from Umit ~ 5.95t to 
U = 6.2t. Second, by examining the high-temperature 
precursors of the transition, we can infer about the ex- 
istence of the first-order transition at large values of U. 
A thermodynamic indicator that signals the proximity to 
the second-order transition line is a peak in the compress- 
ibility {dn/dfi)\T as a function of doping that develops 
at temperatures well above Tcr and reaches its maximum 
at Tcr (see also discussion of Fig. [16]). We found this 
precursor of Tcr for all values of C/ > Umit at progres- 
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sively larger doping as U increases. Another signature 
that can be identified as a precursor of the transition 
is the approximate vanishing of the expansion coefficient 
{dn/dT)\^ = (crossing point of the isotherms in Fig. 2]). 
This signal occurs close to the spinodal line fic2{T) and 
extends to temperature well above Tcr. This feature oc- 
curs for all values of J7 > Uc2 and at progressively larger 
doping as U increases. 



Our results show these high temperature precursory 
signatures of the transition (other indicators will be dis- 
cussed is Sec. IVIII[) for all values of the correlation 
strength U > Umit and hence point toward a critical 
temperature decreasing with U and occurring at progres- 
sively larger doping as U increases. We estimate that the 
critical doping at which Tci- occurs moves from zero dop- 
ing at the Mott point Umit ~ 5.95i, to n ~ 9.95 at 
U = 6.2t, and to w 0.88 at U ^ 12. Ot. Some of the ther- 
modynamic effects close to Tcr will be analyzed in more 
detail in the next section. Here we note that the physics 
arising from the second-order critical line strongly mod- 
ifies not only the thermodynamic response, but also the 
dynamics: as we shall see in Sec. I VII there is a funnel- 
shaped region above the critical line where scattering is 
anomalously large. 



Such precursory signals strongly indicate the occur- 
rence of the first-order transition surface ending at a finite 
temperature second-order transition line for U > Umit 
and thus also for U > 6.2t for temperature lower than the 
range we explored. In addition, these precursory signals 
take place not only for all values oi U > Umit, but also 
without apparent qualitative change. This suggests to 
disfavor other interpretations where a qualitative change 
in these high temperature signatures is expected, like the 
possibility that the critical line ending our first order 
transition becomes a quantum critical line or point^liZ^, 
or a tricritical point at some U. 



Additional support for the approach of a finite T 
second-order transition line follows from the critical slow- 
ing down observed in our CTQMC calculations close to 
that line. Critical slowing down is a widespread and 
standard indicator— that the system is near a critical 
threshold. We find enhancement of Monte Carlo fluctu- 
ations and also enhancement of the number of the cel- 
lular DMFT iterations to attain self-consistency as the 
model parameters approach Tcr- To obtain reliable re- 
sults, a large number of Monte Carlo samplings and of 
the order of hundreds of iterations are used in our nu- 
merical simulations. This effect is well documented in 
single-site DMFT studies of the Mott transition^iiil. In 
cellular DMFT investigations, it is reported in Ref. [s^ 
at i7 = 12. around n ~ 0.88; here we detect this phe- 
nomenon down to Umit at progressively smaller doping, 
hence linking the critical slowing down to Mott physics. 



V. THERMODYNAMICS OF THE 
TRANSITION 

Having established the first-order character of the tran- 
sition, we now focus on the thermodynamic properties of 
the phases that this transition separates. We discuss four 
basic features. 

First we prove that these phases are thermodynami- 
cally stable. We start from the fundamental thermody- 
namic relation 



dE{S, n, U) = TdS + ndn + DdU, 



(7) 



where it is understood that energy entropy S*, number 
of particle n and double occupancy D are normalized 
per lattice site. The model parameters are T, /i, [/, so it 
is natural to use the grand canonical potential $7. The 
appropriate Legendre transform leads to 

dVl = d{E -TS - nn) = -SdT - ndfi + DdU. (8) 

From the expression for dil in Eq. ([5]) and dU = from 
now on, we have 



on 

df 



(9) 



The concavity of the entropy implies that the grand 
canonical potential is minimum at equilibrium. Thus: 



^ [dT d^i) 



{dTf 



\dTdfi 




(dT dfi) < 0. 



(10) 



If we consider a constant T plane the last relation is 
obeyed if 



(11) 



i.e. if the compressibility is positive. In the general case, 
we need to add the requirement that the eigenvalues of 
the above matrix be negative, i.e that the determinant 
be positive: 



\dTdfiJ 



> 0. 



Notice that the specific heat 



T 



dS_ 
df 



(12) 



(13) 



must be positive as a consequence of the above two in- 
equalities. Given 



/ d^n\ _ /dn 
\dTdfi) ^ ^ \df 



(14) 
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n 


/ a \ 
I — 1 

ydfj, J J, 




iar) fx 


det 


0.982 


0.201 


-0.632 


-1.237 


4.420 


0.971 


0.272 


-0.943 


-1.899 


6.944 


0.941 


1.146 


0.116 


0.585 


3.170 


0.921 


0.128 


0.055 


0.648 


4.725 



TABLE I: Particle density n, compressibility dn/dfi\T, expan- 
sion coefficient dn/dT\^, derivative of total energy dE/dT\fi 
and the determinant Eq. 1121 The first (last) two lines cor- 
respond to values of n in the UD (OD) phase. Other data 
parameters are: U = 6.2t > Umit, and the derivatives are 
taken at T/t = 1/70 which lies below Ta- The determinant 
is positive, implying the thermodynamic stability of both the 
UD and OD phases. 



and the expression for dE, Eq. [71 we can extract the 
specific heat from our data since it has the equivalent 
expression: 



T 



dS_ 
df 



dE 
df 



dn 
df 



(15) 



Clearly Fig. 2] shows that dn/dji > 0. As an example, 
we display these quantities for a few points in Table I to 
demonstrate that the stability condition is satisfied. The 
data are for U = 6.2t > Umit and T/t = 1/70 which lies 
below Tcr ~ 1/65. The determinant Eg. 1121 and dn/d^ 
are positive for all values displayed. 

Second, the surprising result that emerges from our 
calculations is the possibility of a first-order transition 
between two paramagnetic metallic phases. The filling 
n is one direct way to discriminate the two phases that 
coexist below the second-order transition line: the OD 
phase at lower filling and the UD phase at higher filling. 

Third, the Clausius-Clapeyron relations reveal that the 
UD phase has smaller entropy and smaller double oc- 
cupancy than the OD phase. The proof works as fol- 
lows. Along the first order transition boundary, the 
grand potential fl is identical in the two phases. This 
means that if we move along the first-order transition 
line on either side, the change in Q is the same in both 
phases: dOuo = d^^OD- If we consider a constant U 
plane [dU — 0) and indicate by subscript c quantities 
measured along the first order transition line, it immedi- 
ately follows that 



dTc _ nun - npn 
dfic Sob — S'ud 



(16) 



Our calculation shows that Tc increases as fic increases 
(i.e. the first-order line bends toward the Mott insulator). 
This implies that the UD phase has a lower entropy than 
the OD phase. In an analogous way, by taking a constant 
T plane, one obtains 



dUc 
dfic 



nuD - npD 
-Dud — -DoD 



(17) 



Our calculations show that fic decreases as Uc increases. 
Hence, the UD phase has lower double occupancy than 



the OD phase. This is as expected and suggests again 
that in the UD phase the correlations are stronger. 

Finally, it is important to stress that the entropy ex- 
hibits a maximum as a function of doping or chemi- 
cal potential in close proximity to the spinodal surface 
fic2{U,T) and its high temperature crossover,. This can 
be seen as follows. Above the critical coupling Uc2{T), 
and for the temperature range explored in our work, the 
isotherms in the n — fj, plane in Fig. |3] approximately 
cross at a finite value of doping. That doping where 
the isotherms cross increases with U. The crossing of 
the isotherms translates into a vanishing expansion co- 
efficient (dn/dT)^ = and thus into an extremum in 
the entropy {dS/dfi)T = since from the grand potential 
Eq. [8] we have the Maxwell relation 



(18) 



As long as the compressibility {dfi/dn)T,u is not sin- 
gular, this also implies {dS/dn)T,u = because from 
S — S{T, /i, U) we can write 

dS = (dS/dT)^ udT + {dS/dfi)T udfi + {dS/dU)T t^dU 

'(19) 

that in turn leads to 

(dS/dn)T,u = {dS/dfi)TM{dfi/dn)T,u- (20) 

The extremum of entropy as a function of filling is re- 
ported also in Refs. 14^16811781 The physical origin of a 
peak of entropy as a function of doping will be discussed 
in Sec. I VIII Here we can anticipate that this feature 
results from reorganization of the basic electronic excita- 
tions of the systems across the OD-UD transition. 



VI. PHASES CHARACTERIZATION: CLUSTER 
QUANTITIES 

Our results show that the normal-state phase diagram 
of the two-dimensional Hubbard model is controlled by a 
first-order transition between a metal and either a Mott 
insulator or another metal that evolves continuously from 
the Mott insulating state. The phases separated by this 
transition have the same symmetry but different thermo- 
dynamic properties, like density, compressibility, entropy, 
double occupancy. 

In this Section we further characterize the nature of the 
sequence of phases that arises upon doping. We discuss 
the cluster Green's function and the cluster self-energy 
that contain direct information about the effects of elec- 
tronic correlations. Our main finding is that critical be- 
havior, as seen in a large scattering rate, originates from 
the spinodal surface Hci{U,T) that delimits the under- 
doped phase. As the carrier concentration moves away 
from that transition, towards both low or high doping, 
coherent electronic behavior is recovered. 
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FIG. 5: Real and imaginary parts of the cluster Green's 
function Giiiiujn) aX U = 7.0t and low temperature T/t — 
1/100 for several dopings. Left panels show the results 
for the orbital K = (0,0), central panels for K = (7r,0) 
and right panels for K — (7r,7r). The finite (zero) value 
of the imaginary part of the cluster Green's function at 
(jj„ — >■ indicates the metallic (insulating) character of the 
solution. Accordingly, three regimes can be distinguished: 
metal (circles), insulator (crosses) and a strong cluster mo- 
mentum differentiation, with orbital K — (vr, tt) gapped and 
K = (0,0), (7r,0) gapless. The first and the last behav- 
ior characterize respectively the overdoped and the under- 
doped phases. The dopings on the figure correspond to 
II = -2.0, -1.5, -1.3, -1.2, -1.15, -1.05, -0.2, 0.0 



A. Green's function 

We begin with a specific example. Fig.[5]shows the real 
and imaginary parts of the cluster Green's function for 
several dopings at U — 7t, above Umit, and at low tem- 
perature T/t = 1/100. As discussed in Sec. HIl in the 2x2 
plaquette scheme the cluster Green's function in cluster 
momentum basis is a diagonal matrix with diagonal en- 



tries G 



(0,0) J G(o,7r), '^{■JT,0), '-f(7r,7r) 



Gr 



G, 



These elements can 



be thought to represent a coarse grained average of the 
lattice Green's function over part of the Brillouin zone. 
The curves display a systematic evolution as the Mott 
insulating state is approached. 

Let us analyze the low frequency behavior of the real 
and imaginary parts of the cluster Green's function. The 
extrapolated value of ImG'K(a; — > 0) is an estimate of 
the local density of states at the Fermi level, a quantity 
accessible in photoemission experiments. For all dopings 
shown, the (tt, 0) component carries most of the spectral 
weight at the Fermi level. At large values of doping, the 
imaginary part is finite in all orbitals (orbitals refer to K 
values), showing that the system is metallic (see curves 
with circles). This behavior distinguishes the overdoped 



phase. On the other hand, at zero doping, the imaginary 
part goes to zero for all orbitals, so it indicates insulating 
Mott behavior for all cluster momenta (see curves with 
crosses). In between these two regions, the imaginary 
part remains finite in the orbital (tt, 0) and (0,0), while 
it shows a clear tendency to go to zero in orbital (tt, tt) 
(see curves with square symbols). This implies that the 
former two orbitals are metallic and the latter is insulat- 
ing, in striking contrast with the OD and non-interacting 
cases. This phenomenon of strong momentum space dif- 
ferentiation^""—!^ characterizes the UD phase and has 
a transparent physical interpretation: doping the Mott 
insulator occurs gradually in certain cluster momenta, 
leading to emergence of metallic behavior in some re- 
gions of the Brillouin zone but not in others. Physi- 
cally, this feature is the coarse-grained manifestation of 
the continuous appearance of the Fermi surface out of 
the Mott insulating state. In the 2x2 plaquette scheme, 
the cluster momentum differentiation in the UD phase 
found in the low frequency behavior of ImG K{'-^n) is 
compatible with arc or pocket formation on the Fermi 
surface. Distinguishing between these scenarios falls out- 
side the scope of the present study. One has to rely on 
periodization schemes to obtain the lattice Green's func- 
tional^""—. Another possibility, recently explored in the 
context of DGA calculations, is to consider progressively 
larger clusters^ii^i^^ . Other band structure effects not 
included here can also be important^. 

Let's now study the real part of the Green's func- 
tion, displayed in the top panels of Fig. [SJ Amongst 
other things, this quantity encodes information about the 
particle-hole symmetry of the orbital at low frequency as 
follows from the spectral representation of the Green's 
function, 



G{iuj„ 



duj ImG^Lj) 



TT lljj„ — OJ 



(21) 



In all the doping range, ReG'(o,o) ('^ 0) is positive, im- 
plying that the largest spectral weight at low frequency 
lies below the Fermi level. The opposite sign is found for 
the orbital (tt, tt), meaning that it has most of the low fre- 
quency spectral weight above the Fermi energy. In sharp 
contrast, the plot of ReG'(7r,o) («w„) reveals a change of 
sign at a characteristic doping. This means particle-hole 
symmetry of the (tt, 0) orbital at that doping, which is re- 
markably close to the doping where the above mentioned 
strong momentum space differentiation in ImG(7r,o) (^'^n) 
sets in. We shall see in Sec. I Villi that this property is 
closely associated with the large scattering rate found 
near the critical transition line. 

So far we have given a qualitative overview of our data 
for a single value of U and a single low temperature that 
is nevertheless higher than first-order transition line Tcr- 
This set of results points to the following picture. Metal- 
lic behavior in all orbitals characterizes overdoped phase. 
Strong modulation of spectral weight at the Fermi level, 
where (7r,0) and (0,0) orbitals are metallic and (tt, tt) is 
gapped or damped, signals the onset of the underdoped 
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FIG. 6: Extrapolated zero frequency value of the imaginary part of the cluster Green's function, -ImGKi^ — 0), as a function 
of doping S = 1~ n, for several values of the interaction strength U. For each figure, the upper panel shows cluster momentum 
K = (7r,0) (circles) and the lower panel displays K = {"KjIt) (triangles). Note the difference of scale of the y-axis. The data 
shown are for temperatures T/t = 1/10 (blue dotted), 1/25 (green dashed), 1/50 (red solid), 1/100 (black dot-dashed). When 
two solutions are found to coexist, the solutions obtained following the metallic and the insulating solution are indicated as full 
and open symbols respectively. This observable measures the density of states at the Fermi energy averaged in a coarse-grained 
cluster momentum region. Doping of the Mott state occurs gradually in cluster momentum, with the metallization that starts 
first in the orbital K — (tt, 0) and (0, 0) while K = (tt, tt) remains insulating. This behavior characterizes the underdoped 
phase and is highlighted by the gray (orange) background for the lowest temperature displayed. The transition between the 
underdoped and overdoped phase can be first order or it can be a crossover depending on U and T. In the latter case, for 
concreteness we define the boundary, illustrated in the figure for T/t = 1/50, by -ImG(^,^)(cj — >■ 0) = 0.05. 



phase. In between these two phases, the (tt, 0) orbital 
becomes strikingly symmetric at low frequency. 

To understand if these are genuine properties of the 
phase diagram, we compile the zero-frequency extrapola- 
tion of the imaginary and real parts of the cluster Green's 
function for a large range of dopings and temperatures. 
We first discuss the behavior of ImGi^ (cj — >■ 0) as a func- 
tion of doping, shown in Figs.|ni[a)-(h) for different values 
of U. To emphasize the cluster momentum differentia- 
tion, each figure has two panels containing data for the 
orbital K = (tt, 0) (curves with circles in the upper panel) 
and (tt, tt) (curves with triangles in the lower panel). The 
(0, 0) orbital behaves similarly to the (tt, 0) one and is not 
shown. The (tt, 0) orbital contains most of the spectral 
weight at the Fermi level and shows the largest temper- 
ature variation. Below the second-order critical line, two 
solutions of the cellular DMFT equations coexist for each 
orbital K. In the overdoped phase ImG/f (w 0) goes to 
a finite value for all orbitals K, revealing the metallic na- 
ture of this phase. On the other hand, in the underdoped 
phase the (tt, 0) and (0, 0) components of ImG/f (w 0) 



go to a sizable value, except asymptotically close to zero 
doping, while the (tt, tt) orbital remains basically gapped 
or with a very small value. The gray (orange) background 
in the lower panels highlights this region for the lowest 
temperature shown. At the spinodal surfaces /id {T, U) 
and fic2iT, U) these two solutions exhibit an abrupt jump 
in each orbitals: from essentially zero to a finite value 
in the case oi K = (7r,7r), and from a finite value to 
a larger value in the case of if = (tt, 0) and (0, 0) (see 
U/t = 5.6,5.8 and 6.0). Above the critical end line the 
two solutions merge and the transition between the two 
phases occurs as a crossover. 

Since Tcr decreases with increasing [/, an immediate 
implication of our results is the following: a study focus- 
ing only on a large interaction strength U will detect the 
changes in ImGi<-(w — t- 0), but not the underlying transi- 
tion between two metals at finite doping, due to the low 
energy scales involved (see U/t = 7.0,9.0,12.0). On the 
contrary, our exhaustive scan of the 3D phase diagram 
allows us to show a direct link between the momentum 
differentiation and the finite doping first-order transition. 
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FIG. 7: ReGif (oj — >■ 0) as a function of doping S, for several values of the interaction strength U. Circles indicate the data for 
the K = (tt, 0) orbital, triangles for the K = (tt, tt). The data shown are for temperatures T/t — 1/10 (blue dotted), 1/25 (green 
dashed), 1/50 (red circles), 1/100 (black dot-dashed). When two solutions are found to coexist, the solutions obtained following 
the metallic and the insulating solution are indicated as full and open symbols respectively. The assignment of symbols is the 
same as in Fig. [6l This observable measures the low frequency asymmetry of the orbitals. The (tt, 0) orbital crosses zero at a 
characteristic doping close to the underdoped-overdoped transition. 



We note that in the UD phase, ImG'(7r.x) ~^ 0) does 
not actually seem to extrapolate exactly to zero, but 
rather to a vanishingly small value. Even if the extrap- 
olated value decreases as T is lowered, and turns from 
sublinear to overlinear dependence on doping, we ratio- 
nalize this feature by the fact that in cellular DMFT the 
orbitals K are coupled through both the Coulomb inter- 
action U and the hopping, making difficult a sharp de- 
coupling among the orbitals at finite temperature. The 
latter behavior is found for example in DCA calcula- 
tion sii i' -SA' ^tSS , where the coupling among the orbitals 
occurs through the interaction U only. 



Fig. [7] shows the extrapolated values of the real part 
of the cluster Green's function, ReGxi^ — s- 0), a mea- 
sure of the particle-hole symmetry of the orbitals at low 
frequency. Here the most striking feature is the change 
of sign in the (tt, 0) component that occurs, apart from 
the obvious particle-hole symmetric line = 0, at 5 = 0, 
at the spinodal surface ^c2iU,T) that signals the disap- 
pearance of the OD phase. This feature is not restricted 
to temperatures below the critical line Tcr but also con- 
tinues for temperatures above it. 



B. Self-energy 

Precious information about the nature of the two 
phases separated by the first-order transition follows from 
the analysis of the cluster self-energy. In Fig.[5]we present 
a selection of self-energies for the same parameter values 
as in Fig. [5j 

An instructive way to look at the data is to consider 
the concept of Fermi liquid coherence within the cluster 
momentum. In Landau's Fermi liquid theory the low en- 
ergy excitations of an interacting system can be described 
as long-lived particles, called quasiparticles. This Fermi 
liquid regime applies below a characteristic temperature, 
the coherence temperature, and a necessary condition for 
its existence is that the imaginary part of the self-energy 
goes to zero as w — ?> with the form S ?s iuj^A + B with 
A negative. 

The self-energy of both the (0,0) and (7r,7r) orbitals 
in Fig. [8] obeys this relation for all the metallic states 
(circles and squares) of the doped Mott insulator, reveal- 
ing Fermi liquid coherence of the electronic excitations. 
Note that the linear low frequency behavior of ImE(o,o) 
and ImE(.n. persists in the Mott insulating state (lines 
with crosses) , where the particles are localized. There the 
system is gapped (see the corresponding low frequency 
vanishing behavior of ImG(o,o)('^n) and ImG(^.7r)(i^n) in 
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FIG. 8: Real and imaginary part of the cluster self-energy 
Sx(iwn) at U = 7.0t and the low temperature T/t = 1/100 
for several dopings. Left panels show the results for the orbital 
K = (0, 0), central panels for K = (tt, 0) and right panels for 
K = (7r,7r). The assignment of symbols is the same as in 
Fig. El 



Fig.O. 

In contrast, the imaginary part of the (tt, 0) self-energy 
shows a non-monotonic behavior as a function of dop- 
ing. The metallic state at both large and small dop- 
ing displays Fermi liquid coherence (ImE(^ o)('^n — > 0) 
goes linearly to zero), but at intermediate doping, at the 
crossover between the underdoped-overdoped phases, the 
large finite intercept indicates a sudden drop of electronic 
coherence (see curves for 6 = 0.05,0.06). Finally, the 
Mott insulating state is reached and a clear pole appears 
in the imaginary part of the (tt, 0) self-energy at 6 — 
and /i = 0. 

The main results of the above analysis are a buildup of 
incoherent metallic behavior close to the first-order tran- 
sition and its high temperature crossover. It is of key 
importance to investigate if this picture survives in the 
whole phase diagram. To this end, we show in Fig. [9] 
the effective chemical potential, which corresponds to 
the chemical potential renormalized by correlation ef- 
fects, and is defined by fi^g — fJ- — ReY^Ki^ — 0). In 
Fig. [TU]we plot the scattering rate F^, estimated from 
the zero-frequency extrapolation of the imaginary part 
of the cluster self-energy = — /mEj(-(a; — )> 0). Although 
these quantities are not strictly equal, we will use them 
interchangeably. We use a linear interpolation and we 
have verified that a quadratic interpolation gives quali- 
tatively similar results. For clarity we show results only 
for the orbitals K = (tt, 0) (circles) and (tt, tt) (triangles). 

The most important result, shown in the upper panels 
of Fig. [TUl is a large scattering rate — /mSj^ o)('^ ~^ 0) 
for a finite range of dopings peaked at the first-order 
transition between two types of metals. A large scat- 



tering rate is in conflict with Fermi liquid theory and 
indicates that short-lived (i.e. incoherent) excitations, 
and not long-lived quasiparticle, produce the metallic be- 
havior in that regime. Nevertheless, if one goes away 
from the transition, the scattering rate rapidly drops 
to small values and thus crosses over to a conventional 
Fermi liquid with small intercept. The overall behavior of 
— /mS(^ o)('^ ~^ 0) unequivocally shows that the source 
of scattering comes from the transition between the un- 
derdoped and the overdoped phase. Some recent studies, 
using different cluster methods, already reported a peak 
in the scattering rate at finite doping, associating it to a 
pseudogap phenomenon^^ , to a competing superexchange 
and Kondo scale^^*^, or to a selective metal-insulator 
transitionSS-jSi.. Our contribution is twofold: first, we link 
this feature to the spinodal 5ci{T) and its high temper- 
ature precursor; second, by tracking this feature in the 
[U —T—5) phase diagram, we show that it is connected to 
the Mott endpoint C/mit, hence unveiling its unexpected 
origin in the Mott physics. 

Let us study in detail the behavior of the self-energy 
in the two phases. The overdoped phase has a rich be- 
havior as doping is reduced towards the first-order tran- 
sition surface. At large doping, far from the transition, 
the effective chemical potential of all orbitals merges, as 
demonstrated in Fig. [9] This is clear at large U (cf. 
Fig.|nUg),(h)) and indeed occurs at smaller values of U for 
similar dopings (not shown). Hence, the self-energy is lo- 
cal (momentum independent) and the single-site DMFT 
solution captures the physics of the problem. As the 
transition is approached, /i^ for different K separate. 
Similar trend has been found in DCA calculations for 
different sizes of cluster~. The corresponding scatter- 
ing rate Tk in Fig. [10] is quite small, decreases with de- 
creasing T and monotonically increases as the doping is 
reduced towards the first-order transition surface. Even- 
tually, at the spinodal 5c2{U, T), E shows a sudden jump 
(Figs. I9ll0p and interestingly the effective chemical po- 
tential of the orbital that carries most of the spectral 
weight at the Fermi level, /i^ff changes sign. The over- 
all low frequency behavior of ImS^ within the overdoped 
phase indicates that the system is a strongly interacting 
Fermi liquid and the effects of electronic correlations dra- 
matically increase as the first-order transition or its high 
temperature crossover are approached. 

Let's now turn to the underdoped phase that origi- 
nates from the Mott insulator. The analysis of the clus- 
ter Green's function revealed that the (tt, tt) orbital is 
gapped but the other orbitals are gapless. Close to the 
Mott insulator and at low enough temperature, the scat- 
tering rate Fj^ o) is small and decreases as T decreases. 
As the the boundary with the first-order transition to 
the overdoped phase is approached, F^^ o) increases and 
eventually reaches its maximum close to the spinodal sur- 
face /ici(J7, T) or near the crossover line above the critical 
line. A further confirmation of this picture emerges from 
the striking temperature evolution of ^(tv,q} displayed in 
Fig. [TUl As the temperature increases, the value of ^{Tr,o) 
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FIG. 9: Effective chemical potential /i^ff = /i — ReT,K{'-o — >■ 0) renormalized by the electronic correlation versus 5. Data are 
shown for the cluster momenta K — {it,0) (circles) and K = {n,n) (triangles) and are obtained for temperatures T/t = 1/10 
(blue dotted line), 1/25 (green dashed line), 1/50 (red solid line), 1/100 (black dot-dashed line). Note that /i'g'^' crosses zero 
at a characteristic doping close to the OD-UD transition, or more precisely at the spinodal I-Lc2{U,T) and its high temperature 
crossover. The hatched grey lines indicate the region where /i^g'^'' exceeds the noninteracting bandwidth of the (tt, tt) orbital. 



at its maximum increases as does its width in doping. 
Interestingly, the lower is Tcr the higher is the tempera- 
ture range at which the value of the (tt, 0) scattering rate 
displays a maximum with doping. Note that the overall 
maximum as a function of U peaks around Uun, relating 
directly this phenomenon to Mott physics. 

In contrast, the (tt, tt) orbital is insulating and there- 
fore r(7r,7r) looses its usual meaning of a scattering rate. 
In search for a possible explanation, we note that, in 

the temperature range explored, At^ff''^^ has a large value, 
that can even exceed the noninteracting bandwidth for 
sufficient large U, as marked by the hatched grey lines 
in Figs. M f,g,li- It has also been discussed, in the con- 
text of DCA framework and for U — 7t, that a pole in 
the self-energy is responsible for the insulating behav- 
ioriliHi69. Our cellular DMFT results for U = U.Ot can 
be compatible with the existence of that pole, however for 
smaller values of U that pole seems not to be present or is 
strongly reduced, as was already pointed out in Ref. [ssl . 
Even if this mechanism remains to be better clarified, 
our contribution is to track the origin of the large scat- 
tering rate to the spinodal surface fici{T, U) and its high 
temperature crossover. 

Finally, we stress that the coherence of the (tt, 0) or- 
bital in the underdoped phase at low doping, as signaled 
by the vanishing of the corresponding scattering rate in 



Fig. [TUl does not imply that the system as a whole is a 
Fermi liquid. As discussed before, the (tt, tt) orbital re- 
mains gapped in this phase suggesting a "small" Fermi 
surface compatible with arc or pocket formation. In other 
words, it is more likely that in the underdoped phase 
electronic coherence does not result from quasiparticle 
propagation but from another mechanism yet to be iden- 
tified. In the next section we shall argue that coherence 
behavior is associated to hole doped carriers moving in a 
background of spin singlets formed by the superexchange 
mechanism. 

What picture emerges from the analysis of cluster 
quantities? Breakdown of Fermi liquid behavior, revealed 
by a dramatic enhancement of scattering rate o)j is 
realized over a large region of the phase diagram close 
to the first-order transition between the underdoped and 
overdoped metals and reaches its overall maximum at the 
Mott endpoint J/mit- Thus our results pin down that in- 
coherent states come from Mott physics, even for regions 
way beyond half filling. As we move away from the tran- 
sition, the lifetime of the low energy excitations of the 
metallic state sharply increases. 

Fermi liquid breakdown beyond half filling is unam- 
biguously tied to Mott physics. The last piece of this 
puzzling scenario is to identify the physical mechanism 
driving this unconventional behavior. This is the focus 
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FIG. 10: Extrapolated zero- frequency value of the imaginary part of the cluster self-energy, -ImE/f (i^ 0) as a function of 
doping 5 and for several values of the interaction U. For each figure, the upper panel shows cluster momentum K = (tt, 0) 
(circles) and the lower panel displays K = (tt, tt) (triangles). This quantity is proportional to the scattering rate Tk- The 
data shown are for temperatures T/t — 1/10 (blue dotted line), 1/25 (green dashed line), 1/50 (red solid line), 1/100 (black 
dot-dashed line). -ImE(^ o)('^ ~^ 0) is peaked close to the spinodal surface fj,ci{U,T) of the first-order transition and to its high 
temperature crossover. That peak reaches its overall maximum in the 3D phase diagram at the Mott end point Umit ~ 5.95 
and S — (see Figs. |10(c)|10(d)[ |. As the temperature increases, the value of -ImE(^ o)('^ ~^ 0) s-t ils maximum increases as it 
does its width in doping. The peak in -ImE(.^^o) (uj — > 0) then disappears at a characteristic temperature which is progressively 
higher as U increases. 



of the next section. 



VII. ROLE OF SHORT-RANGE SPIN 
CORRELATIONS 

To shed light on the physical origin of the incoherent 
metallic behavior associated with the first-order transi- 
tion, we now turn to the analysis of the short-range spin 
correlations. In fact, cellular DMFT takes into account 
on equal footing local quantum fluctuations and finite 
length spatial fluctuations (within the cluster size, here 
a 2 X 2 plaquette). The inclusion of spatial correlations 
manifests itself mathematically in the non locality of the 
self-energyiSi3S. In this section we first compare our cellu- 
lar DMFT phase diagram with that of single-site DMFT, 
where magnetic correlations in space are absent. The dif- 
ferences between the two phase diagrams directly reveals 
the effects of the short-range correlations. Then we char- 
acterize the phases that occur across the transition in 
terms of these fluctuations. As we shall see, in the Mott 
insulating state the electrons are localized due to strong 
screened Coulomb interaction and their spins lock into 
singlets states due to the superexchange mechanism. In 



the UD phase, i.e. the metallic state that evolves out 
of the Mott insulator, the system is still dominated by 
singlet states and metallic behavior is associated with 
doped charge carriers that propagate in a background 
of spins with strong singlet correlations due to superex- 
change. Finally, in the OD phase, the basic excitations of 
the system are quasiparticles due to conventional Kondo 
screening. There are analogies with the suggestion of 
Ref. M. 



A. Comparison with the single-site DMFT phase 
diagram 

A natural way to understand the effects of the short- 
range spin correlations on the Mott transition is to com- 
pare our cellular DMFT phase diagram, Figl^l that 
accounts for these magnetic correlations, with that of 
single-site DMFT, where those correlations are absent. 

Sustained investigations have firmly set the single-site 
scenario for the Mott transition in the single band Hub- 
bard modeli^ii^"— ii^Ti^ and the resulting phase diagram 
as a function of chemical potential, interaction strength 
and temperature can be found for instance in Ref. |24| . 
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Here the key phenomenon that governs the whole phase 
diagram is a finite temperature first-order transition orig- 
inating at the finite temperature Mott critical point. This 
transition occurs between a Mott insulator and a corre- 
lated metallic state, and can take place as a function of 
either U, T or ^. 

Similarly to the single-site DMFT first-order 
transition dominates the normal-state phase diagram also 
in the cellular DMFT case, as we demonstrated in the 
previous sections. But the nature of this first-order tran- 
sition is different and the fundamental difference comes 
from short-range spin correlations considered in cellular 
DMFT. 

Three main effects can be identified. First, they reduce 
the value of the Mott endpoint Umvt- This is as expected 
because the threshold to open a Mott gap increases as the 
frustration at short-distance increase o^^i^^'^"^ . 

Second, the short-range magnetic correlations change 
the shape of the spinodal lines^i^S. In single-site DMFT 
the boundary of the first-order transition bends towards 
the correlated metallic state (as a function of either the 
interaction strength U or the chemical potential /x), so 
that the insulating phase can be reached upon increas- 
ing the temperature of the system. According to the 
Clausius-Clapeyron equation this comes from the fact 
that the correlated metallic state has lower entropy than 
the Mott insulating state. This behavior occurs in several 
correlated systems, like close to the first-order transition 
line in V20j^^^. Physically, this is because in taking 
the large lattice connectivity limit to justify single-site 
DMFT, the magnetic exchange coupling between two 
neighboring spins drops out, thus resulting in a para- 
magnetic Mott insulator with large spin degeneracy of 
order N log 2 (where N is the number of lattice site)^^. 
As discussed in section IVl the opposite behavior occurs 
in cellular DMFT, where the first-order transition bends 
towards the Mott insulating state. This reveals that the 
UD state (and its parent Mott insulating phase) has lower 
entropy than the OD metallic state as a consequence of 
short-range singlet formation that lifts the 2^ degener- 
acy due to the spins of the localized electrons^^fiS. The 
shape of the phase boundary of cellular DMFT at half- 
filling is observed in low-dimensional correlated systems, 
such as the organic conductorsSMS. 

The third and arguably the most surprising effect 
brought about by the short-range magnetic correlations 
is the mere appearance of the UD phase in the cellular 
DMFT phase diagram. In both single-site and cellular 
DMFT solutions, the Hc2{U, T) surface marks the vanish- 
ing of a correlated metallic solution. On the other hand, 
the spinodal fXciiU, T) always marks the vanishing of the 
insulating solution within single-site DMFT, whereas - 
as we demonstrated in Sec. HVl-, this is generally not 
the case within cellular DMFT. Doping instead proceeds 
gradually in certain cluster momenta. We stress that the 
apparent gradual doping of the Mott insulator observed 
over an extremely narrow range of dopings ( 0.002, i.e. 
at least 10 times smaller than what we find) in single- 



site DMFT is purely a finite temperature effect. In our 
case, short-range magnetic correlations do have the ef- 
fect of shifting the spinodal surface fj,ci{U,T) to finite 
values of doping, implying that the first order transition 
can take place between two metallic phases (the UD and 
OD phases) instead of between an insulator and a metal. 
Thus the UD phase occurs only in the cluster DMFT case 
and traces its roots to the short range-singlet formation. 



B. Plaquette eigenstates 

Having established the basic effects of the short range- 
spatial correlations on the structure of the phase dia- 
gram, we now concentrate in the relative importance of 
various many body states in the cellular DMFT solution 
of the Hubbard model. To address this issue we scruti- 
nize the relative statistical weight of the eigenstates |m) 
of the plaquette Hamiltonian i/piaquctto- In the context 
of DMFT calculations, the analysis of the eigenstates of 
the quantum impurity problem has proven to lead to use- 
ful insights on the role of short-range spatial magnetic 
correlation a^^i'^^'^^i^-'^'^^i^^ and also helped to construct a 
physical picture of the nature of the doped carrier o^^'^^ . 

Cellular DMFT maps the lattice onto a 2 x 2 cluster 
of sites immersed in a self-consistent bath. The quan- 
tum impurity problem defined by Eq. ([2|) has a simple 
physical interpretation: as a function of imaginary time, 
the plaquette undergoes transitions between the 4^^ possi- 
ble quantum states by exchanging electrons with the rest 
of the lattice represented by the surrounding bath. The 
bath hybridization function A(w) encodes the dynamics 
of these processes. It is possible to compute the statis- 
tical weight Pfn of the eigenstates |m) of the plaquette 
(see Refs. [6l| for a rigorous definition). It can be inter- 
preted as the relative time that the plaquette spends in 
the cluster eigenstate \m) and is the quantity we focus 
on. 

The plaquette eigenstates \m) can be labeled by their 
number of electrons N, their total spin S and their clus- 
ter momentum K and physically can be thought of as 
representing a coarse grained approximation of the many 
body excitations of the system. 

In Fig. [11] we plot the probability of the plaquette 
eigenstates for U ~ 12. Ot for values representative of 
the Mott insulating state (Fig. [TlT a)). the UD phase 
(Fig.lTlTb)). and the OD phase (Fig. [TIKc)). Results are 
shown for the high temperature T/t = 1/10 and the low 
temperature T/t = 1/50 (left and right respectively). Fo- 
cusing on the region of parameters close to the first-order 
transition, we find that there are only a few plaquette 
eigenstates with large probability. They are the singlet 
with four electrons in the cluster momentum K = (0, 0) 
(black portions of the diagram): 

lAT = 4,5 = 0,ii: = (0,0)), (22) 

the triplet with four electrons and K — (tt, tt) (red por- 
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a) Mott insulating pliase 




b) Underdoped pliase 




c) Overdoped pliase 
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FIG. 11: Statistical weight P of the foUowing cluster eigen- 
states: the singlet S \N = A; S ^ K = (0,0)), triplet T 
\N = A;S = 1;K = (7r,7r)), and doublet D |iV = 3; S = 
1/2; K = (7r,0)), where N, S, and K are the number of elec- 
trons, the total spin, and the cluster momentum of the cluster 
eigenstate. The remnant statistical weight is summed up in 
the histogram R. Data are obtained at [/ = 12. Of and shown 
at T/1 = 1/10 (left) and 1/50 (right). The value of doping 
is (a) 5 = 0, representative of the Mott insulating phase, (b) 
S = 0.02, corresponding to the UD phase and (c) 5 = 0.25, in 
the OD phase. 

tion of the diagrams): 

\N = i,S=l,K={TT,7r)), (23) 

and the doublet with three electrons and K = (tt, 0) and 
the degenerate K = (0, tt) (blue portion of the diagrams). 

|iV = 3,5 = l/2,is: = (7r,0)). (24) 

We sum up the probability of the remnant states (green 
portion of the diagrams) . 

The distribution of the statistical weight in the Mott 
insulator and UD phase is similar as does its evolution 
with temperature. This is expected since the UD phase 
evolves out of the mother Mott insulating phase. Let us 
now summarize the three main differences between the 
two metallic phases. First, in the UD phase the plaque- 
tte electrons are locked into one prevailing configuration 
(half- filled singlet), whereas in the OD phase they fluctu- 
ate among several states. Second, the temperature evo- 
lution of the plaquette eigenstates qualitatively differs in 
the two phases. While charge excitations are weakly T 
dependent in both phases, the basic effect of decreasing 



the temperature strongly suppresses the spin triplet exci- 
tations in favor of the spin singlet ones, especially in the 
UD phase. This has to be associated with the T depen- 
dence of the local density of states at Fermi level shown 
in Fig.[6l Third, the mechanism behind coherence is also 
qualitatively different in the two phases. In the OD phase 
we associate the onset of coherence to standard Kondo 
screening of the conduction electrons of the bath in which 
the plaquette is immersed. In contrast, coherent propa- 
gation in the UD phase occurs when doped holes move 
in a magnetically well defined environment characterized 
by spins bounded into short-range singlet correlations via 
the antiferromagnetic superexchange mechanism. 

Taking the perspective of the quantum impurity model 
immersed in a self-consistent bath is the best way to show 
this point. In Fig. [T^] we show the statistical weight of 
the main plaquette states as a function of doping and 
for several temperatures. In Fig. [T3] we show the ex- 
trapolated zero-frequency value of the imaginary part of 
the cluster hybridization function as a function of doping 
and for different temperatures. The onset of the Mott 
insulating state (5 = for [/ > Umit) is signaled by 
a dramatic growth of the half filled singlet probability 
(black circles). In the Mott state the half- filled singlet 
dominates and its probability increases with increasing 
interaction strength U. The cluster hybridization func- 
tion represents the effect of the bath in which the plaque- 
tte is immersed and with which electrons are exchanged. 
It is gapped at the Fermi level for all cluster momenta 
(i.e. ImAi<-(w — 0) 0). Together, the insulating char- 
acter of the reservoir and the single dominant plaquette 
eigenstate indicate that the electrons of the lattice are 
localized and bound into short-range singlets due to su- 
perexchange interaction. 

In the UD phase, the 2x2 plaquette is still mostly 
occupied by the half-filled singlet (black squares), reveal- 
ing that this phase evolves out of the parent Mott insu- 
lator. However, doping the Mott state populates other 
plaquette eigenstates and the system begins to fiuctuate 
between the half-filled singlet and other configurations. 
Not surprisingly, the primary fiuctuations introduced by 
doping are charge excitations, as shown by the increase of 
the weight of the doublet at = 3 (blue circle) with in- 
creasing doping. In a smaller measure, the probability of 
the spin triplet excitations (red triangles) also grows with 
doping. Interestingly, spin excitations (both of singlet 
and triplet type) show a larger temperature dependence 
than the charge excitations. Note also that the proba- 
bility of the half-filled singlet saturates below a charac- 
teristic temperature. The previous sections have demon- 
strated the metallic, compressible character of the UD 
phase. The fluctuations among various plaquette config- 
urations are allowed by the self-consistent bath, which in 
this phase is metallic (Fig. [T5|) and depends only weakly 
on the cluster momentum, consistent with singlets states 
within the plaquette that dominate the behavior of the 
system. 

Our results thus suggest to consider the UD phase as 
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FIG. 12: Statistical weight P of the following cluster eigenstates as a function of doping 5, for several values of U : the singlet 
|iV = 4,5" = 0,_ft: = (0,0)), the triplet \N ^ 4,S ^ 1,K = (7r,7r)), and the doublet \N = 3,S ^ 1/2, K = (tt.O)) (black 
squares, red triangles, and blue circles respectively) where A'^, S and K are the number of electrons, the total spin and the 
cluster momentum of the cluster eigenstate. The data shown are for temperatures T/t = 1/10, 1/25, 1/50 (dotted, dashed and 
solid line respectively). The statistical weight P of the half-filled singlet crosses 1/2 (green solid horizontal line) roughly at the 
transition between the UD and OD phase, suggesting competing singlet spin versus other spins and charge fluctuations at the 
transition. 



a phase in which the doped-hole carriers propagate into 
a sea of spins bounded into a = 4 short-range singlet 
state (a similar interpretation for low do ping phase of 
the t — J model was suggested in Refs. 53887^. We can 
therefore rationalize why this phase is coherent at low T, 
as revealed by the small scattering rate in Fig. (TU] as 
soon as the spins pair up into singlet pairs, the doped 
holes have no difficulty to delocalize coherently into the 
background of singlet states. 

We now examine the OD phase. Here there is no 
dominant plaquette eigenstate, and the electrons within 
the plaquette fluctuate among several states. Electrons 
spend comparable time visiting the doublet at iV = 4, 
the half- filled singlet and triplet. Also important are 
other types of excitations. At large doping, the charge 
excitations becomes important, so the probability of the 
doublet at = 3 is the largest. As the transition is 
approached, the weight of these excitations decreases, 
while that of the singlet increases. The temperature de- 
pendence of the plaquette eigenstates is weak and mostly 
confined to excitations of spins (both singlet and triplet). 
To complete our survey of the OD phase we analyze the 
effect of the bath in which the plaquette is immersed. 
Fig. [T21 We find that the cluster hybridization function 



displays metallic behavior and weak cluster momentum 
differentiation, especially as the finite-doping transition 
is approached. These results for the OD phase support 
a qualitatively different picture compared with the UD 
phase. At the level of the quantum cluster impurity prob- 
lem, the OD phase is characterized by the electrons of the 
plaquette fluctuating among several conflgurations in a 
conducting bath. Here the bath plays a central role, since 
the electrons of the metallic bath screens the fluctuating 
moments via the Kondo effect. The objects that delo- 
calize in the lattice are both spin and charge excitations 
compatible with Landau quasiparticles that carry both 
spin and charge. As a further support to this picture, we 
stress that at large doping, far from the transition, we re- 
cover the results of single-site DMFT, where the spatial 
fluctuations are absent and where the Kondo screening 
mechanism controls the physics of the model. 

The central question then is: how can the system go 
from the OD phase to the UD phase? At the level of the 
plaquette, the transition between the two phases results 
from the competition between singlet spin excitations 
and spin triplet plus charge excitations. A simple way to 
quantify the increased importance of spin triplet-charge 
fluctuations at the expense of the singlet across the tran- 
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FIG. 13: Extrapolated zero-frequency value of the imaginary part of the cluster hybridization function, -ImAif(cj — >■ 0), for 
cluster momentum K — (tt, 0) (circles) and K — (tt, tt) (triangles) as a function of doping, for different values of the interaction 
strength U. The data shown are for temperatures T/t = 1/10 (blue dotted line) 1/25 (green dashed line), 1/50 (full red line), 
1/100 (black dot-dashed line). This observable measures the density of states of the conduction bath at the Fermi level averaged 
in the coarse grained cluster momentum region. Note the weakly cluster momentum dependence in the underdoped phase, to 
be contrasted with the strong cluster momentum differentiation of both the cluster Green's function and self-energy. 



sition, is to monitor when the probability of the singlet 
state equals the probability of other plaquette states, i.e. 
when Pg = 1/2 (green solid line in Fig. [T^ . This occurs 
roughly at the UD-OD transition. 

The competing spin-charge excitation close to the tran- 
sition was already reported in our recent study, Ref. [tO. 
There, we emphasized the interplay between the doublet 
and the spin excitations. Close to the critical coupling 
C^MiT (cf. Fig.fT^.c.d.e) we linked the coexistence region 
to the crossing of the probability between the triplet and 
the doublet states. In the present study, which explores 
a larger region of interaction strength U and temper- 
atures, we rather highlight the competing singlet spin 
versus triplet-spin and charge fluctuations. This aims to 
refine our previous study because the crossing Ps ~ 1/2 
tracks both the first-order transition at low temperatures 
and its high temperature crossover providing a more gen- 
eral framework for interpreting our results in the whole 
range of U and T explored. 

We can now address the crucial question of the origin 
of large scattering rate close the the /ici {T, U) surface 
and its high temperature crossover analyzed in Sec. I VII 
The results of the plaquette eigenstates describe compet- 
ing^ superexchange and Kondo interactions that give 
rise to unquenched short-range magnetic fluctuations of 
the electrons within the plaquette. This identifies the 
source of scattering close to the transition between the 



UD and OD phase. 

The above description takes the point of view of the 
cluster coupled to a bath. The cellular DMFT results 
have also a complementary interpretation in terms of the 
lattice. In the OD phase the electrons on the lattice are 
coherent because they fill the available k states as in a 
Fermi liquid^-. At zero temperature, this quenches the 
entropy. Here the k space picture is appropriate. On the 
other hand, in the UD phase the electrons pair up into 
short-range singlets and doped charge carriers can move 
coherently into a background of singlet bonds. Here a 
real-space picture is more appropriate. Starting at large 
doping and reducing S towards the transition to the UD 
phase, the electrons reorganize their motion according 
to the strength of short-range magnetic spin correlation 
that increase as the Mott state at (5 = is approached. 
Breaking of the electrons as entities carrying both spin 
and charge occurs at the finite doping transition. This 
releases entropy as confirmed in Sec. |Vl Once the new 
excitations in the form of short range singlets are formed, 
and the system get settled into a magnetically stable con- 
figuration at low T, the system can acquire again coher- 
ence. This perspective identifies the reorganization of 
the basic excitations of the system at the transition be- 
tween the OD and UD phase as the source of scattering. 
Such reorganization also occurs in some recent theoreti- 
cal proposals that envision a fractionalisation of the elec- 
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trons22,i22,. A different type of reorganization occurs in 
the charge 2e boson theorySI. 

Further support for our scenario follows from 
the the zero frequency cluster susceptibility xo = 
/p^(S'^(r)S'J(0))dr, where is the projection of the total 
spin along the z direction of the 2x2 plaquette, shown 
in Fig. [U as a function of doping for several values of 
U and T. In the Mott insulating phase this quantity 
vanishes as temperature is lowered, revealing that the 
system is locked into 5 = singlet states due to the su- 
perexchange interaction. This is in sharp contrast with 
single-site DMFT where the corresponding local spin sus- 
ceptibility diverges in the Mott state. As the doping is 
increased, xo grows approximate linearly with doping. 
This characterizes the UD phase and reinforces the view 
of this region as a phase with hole-doped carriers prop- 
agating in a background of singlet states. On the other 
hand, starting at large doping, a decrease of 5 towards 
the transition has again the effect of slightly increas- 
ing Xo- The mechanism however is different: spin and 
charge excitations are Kondo screened by the conduction 
electrons of the bath. In a two-dimensional Fermi liq- 
uid with spherical Fermi surface, the spin susceptibility 
depends only on the density of states and is indepen- 
dent of doping in the non-interacting limit. Here, as the 
transition is approached with decreasing doping, xo in- 
creases slightly, then approximately saturates or slightly 
decreases, revealing the increased importance of the un- 
quenched spin fluctuations in the plaquette as the tran- 
sition is approached. Notice also the opposite temper- 
ature behavior of xo within the two phases, again sup- 
porting the two different low lying elementary excitations 
in these two regions. In the UD phase in particular, the 
spin susceptibility decreases rapidly with temperature, 
in analogy with the spin-gap phenomenon found in high- 
temperature superconductors2ii22. 



VIII. DISCUSSION 

Understanding the origin of the unusual behavior of 
copper oxide superconductors is a fascinating challenge 
for theorists and experimentalists alike. Recent exper- 
iments are focusing on the normal state phase of these 
systems and are providing new clues on the critical be- 
havior observed between the metallic state at high dop- 
ing and the Mott insulator at zero doping^. In an at- 
tempt to obtain a unified picture of these systems, we 
have mapped out the normal state phase diagram of the 
2D Hubbard model, which is arguably the simplest model 
able to capture the essential physics of doped Mott in- 
sulator^. Cellular DMFT on a 2 x 2 plaquette is our 
theoretical tool to elucidate the phase diagram, allowing 
us to consider both temporal and short-range spatial fluc- 
tuations. Our results provide a unified low-temperature 
landscape where, at the lowest temperatures, the system 
evolves, as a function of the carrier concentration, from 
a Mott insulator at zero doping, through a coherent elec- 



tron liquid in the UD region, and, with further doping, 
to a conventional metallic state in the OD region. Our 
key finding is the first-order transition that takes place 
at low temperature between these two metallic phases at 
finite doping, which is a source of large scattering rate at 
the temperatures relevant for experiment. 

In this section we further discuss the emerging phase 
diagram with a focus on the typical signatures and pre- 
cursors of the first-order transition. By tracking these 
indicators down to the Mott critical endpoint Umit we 
reveal that Mott physics is at the origin of this transition. 
Subsequently we conjecture that the large scattering aris- 
ing from this transition may provide an alternative to the 
quantum critical behavior surmised in the copper oxide 
superconductors22i2i. 



A. Mott physics beyond half-filling 

Fig. [15] shows the temperature-doping phase diagram 
for U = 6.2t representative of the regime U > C/mit and 
summarizes our results for the doping induced transitions 
starting from the Mott insulator. In the T — S plane there 
are five regions of interest: Mott insulator, underdoped 
phase, hysteretic region, overdoped phase, and the non- 
Fermi liquid region arising from the finite-temperature 
crossover between the UD and OD phases. 

Mott insulator. At zero doping the system is a Mott 
insulator, as revealed by a plateau in the occupation at 
the (odd) number n = 1 as a function of fj, (see Sec. IIII[) . 
Here the spins are bound into short-range singlets due to 
superexchange mechanism (see Sec. IVIip . 

Underdoped phase. Doping of the MI state proceeds 
gradually in certain cluster orbitals, with carriers going 
into the (0,0), (0,7r) and the degenerate (7r,0) orbital, 
but not in the (tt, tt) one, that remains insulating. This 
strong orbital differentiation is robust with respect to an 
increase of the size of the cluster immersed in the self- 
consistent mediuni^^. Notice however that in CTQMC 
calculations this gradual doping in cluster orbitals is as- 
sociated with a large orbital differentiation of the cluster 
self-energy, while exact diagonalization methods show the 
latter phenomenon without sign of the former, at least 
for a 2 X 2 plaquette^i^. The basic excitations of this 
phase show coherent behavior, as indicated by the linear 
low frequency behavior of ImE(7r.o) and ImS(o.o)- How- 
ever the system falls outside Fermi liquid theory since the 
(tt, tt) orbital remains gapped in this phase suggesting a 
"small" Fermi surface compatible with the presence of 
arcs or pockets (see Sec. IVip . Furthermore, the objects 
that delocalize through the lattice are holes propagating 
into a background of short-range singlet states, as sug- 
gested by the analysis of the plaquette eigenstates (see 
Sec. lVII[) . The metallic character of this phase is revealed 
by the finite spectral weight at the Fermi level in the 
(tt, 0) and (0, 0) orbitals and by the finite charge com- 
pressibility deduced from the dn/d^lx behavior. Since 
the basic excitations of this phase are not the usual Lan- 
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FIG. 14: Local spin susceptibility xo versus 5 for several values of U. The data shown are for temperatures T/t = 
1/10, 1/25, 1/50 (blue dotted, green dashed, red solid line respectively). In the underdoped phase xo is approximately linear 
with doping, as highlighted by the solid cyan line, substantiating the idea that the doped carriers propagate among short-range 
singlets formed by the superexchange interaction. 



dau quasiparticles, it will be interesting to explore the 
way these excitations carry heat and charge. 

Overdoped phase. Upon further doping, there is a first- 
order transition to a phase that appears to be a conven- 
tional Fermi liquid phase. 

Coexistence region. Below the characteristic tempera- 
ture Tcr, there is an hysteretic region in the T — ^ plane 
where the OD and UD phases are both present. There- 
fore within this region a first-order transition takes place 
where the free energy of the UD and OD phase crosses. 
One of the two phases is metastable, but a computation 
of the free energy to tell which one is metastable would be 
prohibitive. In the T — 5 plane there is forbidden region 
delimited by the spinodals, which corresponds to the dis- 
continuous jump of the occupation n(/i). If we were able 
to exclude the metastable phases, this forbidden region 
would be wider and delimited by the so-called binodals. 

Crossover region. Above Tcr, between the UD and the 
OD phase there is an incoherent phase signalled by a 
large scattering rate in the (tt, 0) orbital. This region has 
a funnel-like shape in the T — 5 parameter space. It is 
noteworthy that the lower is Tcr, the higher is its exten- 
sion in temperature. Future work should better clarify 
the thermodynamic critical behavior of this region as well 
as its signatures on transport properties. We note never- 
theless that there is a maximum in entropy as a function 
of doping close to the boundary of the OD phase, see.lVl 



This occurs because the two phases have different elec- 
tronic excitations. At the first order transition, or just 
above, the electronic excitations of the system are thus 
reorganized leading to a large entropy. 

Let us now review the phase boundaries. As it ap- 
pears from Fig. [131 the phase diagram is dominated by 
the critical point at the characteristic temperature Tcr 
(star symbol) . Below this critical point the metallic state 
separates into two distinct metallic states: the UD state 
to the left of the critical point and the OD state to the 
right. The spinodal lines 5ci and (5c2 emanate from the 
critical point (solid and dashed blue curves with circles) 
and are defined as the loci in the phase diagram where 
the occupation n = \ — 5 undergoes a sudden drop and 
an upward jump respectively. Hence, the region between 
(5c 1 and 5c2 is forbidden and translates into a coexistence 
of phases with different densities in the T — jj, plane with 
hysteretic behavior clearly apparent in the n(/Lt) curves 
(see Sec. lIVp . Above the critical point at Tcr there is a 
characteristic crossover line, denoted by the dotted blue 
line with circles, and defined by the peak in the charge 
compressibility dn/dfj,\T. 

To detect directly the first-order transition by looking 
for hysteresis in some observables, like the occupation, 
requires access to low temperatures and a careful scan 
of phase space. Since rapidly moves to low T as U 
increases it is very helpful to find some precursory signals 
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that allow to infer such a transition. 

The first indicators that the system is approaching Tcr 
is the above-mentioned crossover line defined by the max- 
imum in the compressibility and the associated critical 
slowing down in our CTQMC calculations close to that 
line (see Sec. IIVI) . Interestingly, we also found precur- 
sors of the spinodal lines Sd and 6c2, where the UD and 
OD phases cease, respectively, to exist. Signs that the 
model is approaching the spinodal Sc2 include the vanish- 
ing of ReG'(^ o)('^ ~^ 0) and of the renormalized chemical 

potential ^[,9°^ and the crossing of the isotherms n(fi) 
(solid, dashed and dot-dashed red curves with squares 
in Fig. [T5|) . These indicators originate at low tempera- 
tures close to the Sc2 hue and persist for temperatures 
higher than Tcr- The first two phenomena reflect a low 
frequency particle-hole symmetry of the orbital (tt, 0) at 
this critical line. 

We now compile the characteristic precursors of the 
spinodal 5ci, shown as green curves with triangles in 
Fig. 1151 The indicators includes the maximum of the 
scattering rate of the (tt, 0) orbital T — — /mI](^^o)('^ ^ 
0) (solid line), the strong orbital differentiation in the 
density of states at the Fermi level, measured by - 
ImG'i<-(a; — > 0) (dot-dashed curve) and the line where the 
probability of the spin-singlet plaquette state that dom- 
inates at half- filling exceeds 1/2. Below Tcr these phe- 
nomena take place close to the spinodal line ^ci, therefore 
directly signalling the transition. Above Tcr, they may 
serve as indicators that the transition is approaching. 
We note that these precursory signals can be detected at 
rather large temperatures. In addition, we have verified 
that as Tcr decreases with increasing U, these indicators 
can actually be detected at higher and higher tempera- 
tures, a definite advantage when it comes to tracking the 
position of Tcr when it becomes very low. This behavior 
is quite surprising and calls for a better understanding of 
the high T crossover above the transition. It seems that 
when Tcr is very low, it effectively acts as a quantum crit- 
ical point. Quantum criticality influences the physics at 
rather high temperature^^"— . 

We now show that all the above indicators are mani- 
festations of Mott physics. If we follow them as a func- 
tion of the correlation strength U (at a given temper- 
ature T/t = 1/50), we obtain the results shown in the 
6 — U plot of Fig. [ini The critical line Tcr (blue curve 
with circles) originates at the Mott endpoint C/mit with 
coordinates {U/t w 5.95,5 = 0) and moves to progres- 
sively larger doping as U is increased, as discussed in 
detail in Sec. IIVI The line tracking the U evolution 
of the maximum in scattering rate of the (tt, 0) orbital 
(green curve with triangles) originates at the spinodal 
Uci at zero doping and departs from zero doping with in- 
creasing U. Due to particle-hole symmetry, the indicator 
ReG(^^o)('^ — > 0) = (red line with squares) is present 
at zero doping for all values of U but, in addition, at the 
spinodal Uc2 the curve bifurcates and a new branch splits 
away from zero doping as U increases. 

The fact that all signatures of the finite doping first- 
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FIG. 15: Phase diagram temperature T versus doping 5 for 
U = 6.2t, representative of the regime U > Umit- The star 
symbols denotes the second-order finite temperature critical 
point, with coordinates (5, T) ~ (0.05, 1/65). Blue lines with 
circles identify the spinodals Sd (full line) and 5c2 (dashed 
line) where the UD and OD phases respectively disappear. 
Dotted blue line with circles denotes the high-temperature 
crossover of the transition measured by a maximum in the 
compressibility {dn/dfijr- The red lines with squares indi- 
cate the vanishing of ReG(^.o)(i^ — >■ 0) (solid line) and of 
/i^fl'"' (dashed line) and the crossing of the isotherms n{fj,) 
(dot-dashed line) and they are early warning indicators of 
the spinodal 5c2{T). The green lines with triangles are pre- 
cursors of the spinodal 5ci{T). Solid line marks the max- 
imum of the scattering rate r(^ o) = —I'n^'^(TT,o){'^ ^ 0). 
The strong momentum differentiation found in ImGis:(aj — > 0) 
disappears at the dot-dashed line. When that differentia- 
tion occurs as a crossover, we defined the boundary by - 
ImG(^_^)(a; — >■ 0) = 0.05. The dot-dashed line corresponds to 
the doping below which the probability of the half filled sin- 
glet state becomes larger than 1/2. Below the critical point 
the metallic state separates into two different metallic states, 
the UD and the OD phase, with a hysteretic change in the 
occupation (the jumps in the occupation define hatched re- 
gion). The lines above the critical point are crossover lines. 
The Mott insulator (MI) exists only on the 5 = line. 



order transition can be traced back to the Mott critical 
point reveals that the first-order transition between the 
UD and the OD phase is associated with Mott physics. 
This in turn implies that the signature of the Mott tran- 
sition in the two-dimensional Hubbard model, normally 
expected to occur in close proximity to half-filling, ex- 
tends way beyond half filling. 

This surprising conclusiort^S emerges naturally from 
the detailed mapping of the 3D normal state phase dia- 
gram. It is indeed quite difficult to derive such a picture 
from the analysis of the doping driven transition at large 
U only. Studies of this type may detect a critical dop- 
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FIG. 16: Characteristic dopings versus U that signal the UD 
to OD transition, measured at the temperature T/t — 1/50. 
The blue line with circles denotes the maximum of the com- 
pressibility {dn/d^)T- The red full line with squares indicates 
the vanishing of ReG'(^ o)('^ ~^ 0). Green line with trian- 
gles denotes the maximum of the scattering rate Tf^ o) =- 
ImE(^ o)(^ ~^ 0). Data points are displayed for the constant 
temperature T/t = 1/50 and thus they measure the features 
of the three-dimensional phase diagram at that temperature. 
The overall maximum (as a function of U , T and 5) of the 
compressibility follows the critical line Tcr that originates at 
Umit, with coordinates ((7, 5, T) « (5.95, 0, 1/12) indicated by 
the vertical arrow. The overall maximum as a function of U , 
T and 5 of the scattering rate r(^_o) follows the 5ci{U, T) sur- 
face and peaks at UmvT- The zero of the extrapolated value 
of real part of the (7r,0) cluster Green's function at a; = 
bifurcates at the line Uc2{T). 



ing and the strong orbital difFerentiation approaching the 
Mott insulator. It would be hard, though, to find that 
the detected critical doping is indeed a high temperature 
signature of a finite T first-order transition, because of 
the low energy scales involved (i.e. the critical line Ta 
moves rapidly to progressively lower temperature as U 
increases). It would also be hard to make an unambigu- 
ous link between that transition at finite doping and the 
Mott endpoint at zero doping. 

The map of the whole normal-state phase diagram, 
along with our discovery of a first-order transition orig- 
inating at the half-filled Mott critical point and mov- 
ing progressively away from half-filling as U increases, 
should help provide a unified picture of a host of previ- 
ously known results and call for their reinterpretation as 
precursory signals of the transition we have revealed. In 
this respect, and accounting for the differences in meth- 
ods, the pseudogap - Fermi liquid transition reported for 
U larger that the Mott endpoint in Ref. as well as the 
crossover observed at large doping in Refs. HgtHBSI are 
likely to be the same as those observed here. Neverthe- 



less, our analysis suggests that the quantum critical point 
hypothesized in Refs. i93n94i is instead a very low but finite 
T critical point that originates from the influence of Mott 
physics away from half filling. We cannot strictly rule out 
that the critical line that ends our first-order surface be- 
comes a quantum critical line or point at some U (see 
theoretical investigations by Imada in Refs. It^ItBI ). how- 
ever this hypothesis is disfavored in our calculations by 
the fact that the same precursory signs of the transition 
occur for all the values of U > Umit considered, with- 
out apparent qualitative change (see Fig. ^E^i . The same 
argument holds to disfavor the hypothesis of a tricritical 
point at some U in the phase diagram. Several earlier 
studies have also shown that the strong orbital difFeren- 
tiation phenomenon is a robust feature of cluster DMFT 
solutions of the two-dimensional Hubbard model^^. But 
it was unclear what is the origin of the strong momen- 
tum differentiation and of its disappearance with dop- 
ing. This phenomenon finds a natural explanation in our 
phase diagram since it marks the end of the UD phase. 

Even if more investigations should be done to clarify 
both the crossover behavior above Tci- and the zero tem- 
perature limit of our phase diagram, our study reveals 
that a first-order transition inhabits and controls much 
of the physics of the normal-state phase diagram of the 
two-dimensional Hubbard model. 



B. Proposal for critical behavior 

We now turn to the relevance of our results for the 
physics of the copper oxide superconductor o^^'^°°' '-^. 
The difficulties to achieve a unified picture are both ex- 
perimental and theoretical. From the theory side, they 
arise mainly because the phase diagram of the cuprates 
defies three paradigms of condensed matter physics: 
band theory fails to explain the Mott insulator at zero 
doping, Bardeen-Cooper-Schrieffer description does not 
describe the unusual superconductivity and the Landau 
Fermi liquid theory does not capture the metal above 
Tc and the pseudogap phase close to the Mott insulator. 
Therefore, understanding the evolution of these systems 
as a function of doping remains a challenging issue. Since 
the parent compounds are Mott insulators, this calls for 
an understanding of the doping-driven MIT arising from 
doping the Mott insulator Models of this phenomenon 
rely on simple Hamiltonians such as the Hubbard model. 
Here the aim is to ascertain to what extent this mini- 
mal model is able to capture the essence of the complex 
physical behavior of the copper oxide superconductors. 
We mention that in this context the optical lattices of 
cold atoms have opened up new perspectives to test these 



assumption! 
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and determine the phase diagram of 



basic Hamiltonians of condensed matter physicsiSl J£S 

From a methodological viewpoint, dynamical mean- 
field theory approaches allow to follow different phases 
as a function of model parameters. In this paper we 
concentrated on the normal state, by studying the cellu- 
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lar DMFT paramagnetic solution of the two-dimensional 
Hubbard model. Thus the main experimental relevance 
of our work consists in elucidating how the proximity to 
the Mott insulator affects the normal state. 

Our key finding, namely the first-order transition 
between two metals, suggests the conjecture that the 
quantum critical behavior surmised in high-temperature 
superconductors2i21 originates from a constant U cut of 
the finite-doping critical line ending the first-order sur- 
face between two types of metals. That transition occurs 
only when U is larger than the critical value Ud for the 
Mott transition and moves to progressively higher doping 
as U increases, reaching 6 w 0.12 at U = 12t, which cor- 
responds to the doping that maximizes the superconduct- 
ing transition temperature^"^. This suggests that our 
first-order transition lies deep below the superconduct- 
ing dome (cf. Refs. IssIIbsIIstI ). In addition, at large U the 
critical point Tcr is at extremely low temperature and the 
smaller Ta the farther its influence, in the form of large 
scattering rate, extends at finite temperature^^"— . Our 
contribution is to pin down the large scattering rate to 
the infiuence of the Mott phenomenon far away from half 
filling: the first-order transition is the source of anoma- 
lous scattering, which displays a funnel shape with in- 
creasing temperature. Such funnel shape scattering sepa- 
rates underdoped from overdoped phases above the max- 
imum of the superconducting dome in high temperature 
superconductors. Therefore the UD and OD phases, that 
we defined with respect to our metal-metal first order 
transition, should also correspond to the underdoped and 
overdoped regions normally defined in high temperature 
superconductors. Given the differences in methods, the 
crossover that we found between UD and OD phases is 
also the one detected for J7 = 8Hn Refs. lijlesl . 

Further work is planned to make the correspondence 
clearer, but at present at the qualitative level several 
of our results are similar to those found in experiment. 
First, anomalous scattering takes place near a critical 
doping. Contrary to other approaches^^^, the source of 
this scattering is not a quantum critical point between 
two ordered phases, nor the proximity to an ordered 
phase, but a finite temperature crossover at finite doping 
coming from Mott physics at work far from half filling. 
Second, above the critical point, the thermal expansion 
coefficient (dn/dT)^ vanishes and equivalently, the en- 
tropy has an extremum {dS/dy^x = 0^- Using the 
Kelvin formulapiii, this suggests that the thermopower 
changes sign at finite doping''^, consistent experimenl^ii^. 
Finally, as observed in the pseudogap phase^iiS^, there is 
a strong decrease with temperature of the spin suscepti- 
bility in our underdoped phase. 



IX. CONCLUSION 

The complete normal-state phase diagram of the two- 
dimensional Hubbard model as a function of tempera- 
ture, doping and interaction strength reveals a surface of 
first-order transition that ends on a critical line and that 
can separate two different types of metals. That transi- 
tion at finite doping between the underdoped and over- 
doped phase is unexpectedly connected to Mott physics, 
since it originates at the half-filled Mott endpoint and 
moves progressively away from half-filling as the interac- 
tion strength increases. The critical line rapidly moves 
to low temperature as interaction strength is increased 
beyond the critical value for the Mott transition, and it 
leaves finite-temperature signatures similar to those that 
have suggested the existence of a quantum critical point 
near optimal doping in high-temperature superconduc- 
tors. Hence we propose that Mott physics leads to a very 
low temperature critical point separating two metals and 
that this can be the source of anomalous normal state 
behavior near optimal doping in high-temperature su- 
perconductors. Even though this critical point cannot be 
accessed because of intervening long-range order, it con- 
trols the physics at sufficiently high temperature in the 
same way that a zero-temperature Fermi-liquid ground 
state controls the normal state Fermi liquid despite the 
fact that the true ground state may have long-range or- 
der. 

Our proposal is clearly in its infancy and further work, 
already in progress, should reveal the precise nature of 
the criticality along with the relationship between the 
pseudogap phase and the UD region of our phase dia- 
gram. In addition, the study of some ordered states, like 
the antiferromagnetic and superconducting ones, should 
provide useful insight to complete the cellular DMFT pic- 
ture of the Mott transition. Nevertheless, our discovery 
sheds new light on the phase diagram of the 2D Hub- 
bard model and highlights that Mott physics can play a 
pivotal role even far away from half filling. 



ACKNOWLEDGMENTS 

We acknowledge Steve Allen for invaluable technical 
help and Patrick Semon and David Senechal for useful 
suggestions and discussions. This work was partially 
supported by FQRNT, by the Tier I Canada Research 
Chair Program (A.-M. S. T.), and by NSF DMR-0746395 
(K. H.). Computational resources were provided by CFI, 
MELS, the RQCHP, and Compute Canada. 



^ N. Doiron-Leyraud, C. Proust, D. LeBoeuf, J. Levallois, 
J.-B. Bonnemaison, R. Liang, D. A. Bonn, W. N. Hardy, 
and L. Taillefer, Nature 447, 565 (2007). 

^ D. LeBoeuf, N. Doiron-Leyraud, J. Levallois, R. Daou, 



J.-B. Bonnemaison, N. E. Hussey, L. Balicas, B. J. 
Ramshaw, R. Liang, D. A. Bonn, et al., Nature 450, 533 
(2007). 

^ C. Jaudet, D. VignoUes, A. Audouard, J. Levallois, 



26 



D. LcBocuf, N. Doiron-Lcyraud, B. VignoUc, M. Nardonc, 
A. Zitourii, R. Liang, ct al., Phys. Rev. Lett. 100, 187005 
(2008). 

" S. E. Sebastian, N. Harrison, E. Palm, T. Murphy, C. H. 

Mielke, R. Liang, D. A. Bonn, W. N. Hardy, and G. G. 

Lonzarich, Nature 454, 200 (2008). 
^ E. A. Yelland, J. Singleton, C. H. Mielke, N. Harrison, 

F. F. Balakirev, B. Dabrowski, and J. R. Cooper, Phys. 
Rev. Lett. 100, 047003 (2008). 

^ A. F. Bangura, J. D. Fletcher, A. Carrington, J. Levallois, 
M. Nardone, B. Vignolle, P. J. Heard, N. Doiron-Leyraud, 
D. LeBoeuf, L. Taillefer, et al., Phys. Rev. Lett. ICQ, 
047004 (2008). 

^ L. Taillefer, Annual Review of Condensed Matter Physics 

1, 51 (2010). 
* P. W. Anderson, Science 302, 89 (2003). 
^ P. W. Anderson, The theory of Superconductivity in the 

High Tc cuprates (Princeton University Press, Princeton, 

1997). 

^° N. F. Mott, Metal-insulator transitions (Taylor & Francis, 
London, 1974). 

M. Imada, A. Fujimori, and Y. Tokura, Rev. Mod. Phys. 
70, 1039 (1998). 

P. A. Lee, N. Nagaosa, and X.-G. Wen, Rev. Mod. Phys. 
78, 17 (2006). 

A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, 
Rev. Mod. Phys. 68, 13 (1996). 

W. Metzncr and D. VoUhardt, Phys. Rev. Lett. 62, 324 
(1989). 

T. Maier, M. Jarrell, T. Pruschke, and M. H. Hettler, Rev. 
Mod. Phys. 77, 1027 (2005). 
1^ G. Kothar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, 
O. ParcoUot, and C. A. Marianetti, Rev. Mod. Phys. 78, 
865 (2006). 

^"^ A.-M. Tremblay, B. Kyung, and D. Senechal, Low Temp. 
Phys. 32, 424 (2006). 

M. PotthofF, Eur. Phys. J. B (France) 32, 429 (2003). 
^® P. Limelette, A. Georges, D. Jerome, P. Wzietek, P. Met- 

calf, and J. M. Honig, Science 302, 89 (2003). 
^° M. J. Rozenberg, X. Y. Zhang, and G. Kothar, Phys. Rev. 

Lett. 69, 1236 (1992). 

M. J. Rozenberg, G. Kothar, and X. Y. Zhang, Phys. Rev. 
B 49, 10181 (1994). 

M. J. Rozenberg, R. Chitra, and G. Kothar, Phys. Rev. 

Lett. 83, 3498 (1999). 

A. Georges and W. Krauth, Phys. Rev. B 48, 7167 (1993). 

G. Kotliar, S. Murthy, and M. J. Rozenberg, Phys. Rev. 
Lett. 89, 046401 (2002). 

G. Kotliar, E. Lange, and M. J. Rozenberg, Phys. Rev. 
Lett. 84, 5180 (2000). 
^® G. Kotliar, The European Physical Journal B - Condensed 
Matter and Complex Systems 11, 27 (1999), ISSN 1434- 
6028. 

2^ P. Phillips, Rev. Mod. Phys. 82, 1719 (2010). 

S. Lefebvre, P. Wzietek, S. Brown, C. Bourbonnais, 

D. Jerome, C. Meziere, M. Fourmigue, and P. Batail, 

Phys. Rev. Lett. 85, 5420 (2000). 
^® P. Limelette, P. Wzietek, S. Florens, A. Georges, T. A. 

Costi, C. Pasquier, D. Jerome, C. Meziere, and P. Batail, 

Phys. Rev. Lett. 91, 016401 (2003). 
^° F. Kagawa, K. Miyagawa, and K. Kanoda, Nature 436, 

534 (2005). 

F. Kagawa, K. Miyagawa, and K. Kanoda, Nat Phys 5, 
880 (2009). 



T. Ghashi, T. Momoi, H. Tsunetsugu, and N. Kawakami, 
Physical Review Letters 100, 076402 (pages 4) (2008). 
A. Liebsch, H. Ishida, and J. Merino, Phys. Rev. B 79, 
195108 (2009). 

3* M. H. Hettler, A. N. Tahvildar-Zadeh, M. Jarrell, T. Pr- 
uschke, and H. R. Krishnamurthy, Phys. Rev. B 58, 

R7475 (1998). 

A. L Lichtenstein and M. L Katsnelson, Phys. Rev. B 62, 

R9283 (2000). 

^® G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, 
Phys. Rev. Lett. 87, 186401 (2001). 
M. Vekic and S. R. White, Phys. Rev. B 47, 1160 (1993). 
H. Park, K. Haule, and G. Kotliar, Phys. Rev. Lett. 101, 

186403 (2008). 

E. Gull, P. Werner, X. Wang, M. Troyer, and A. J. Millis, 
Europhys. Lett. 84, 37009 (2008). 
'*° M. Balzer, B. Kyung, D. Senechal, A.-M. S. Tremblay, 
and M. Potthoff, Europhys. Lett. 85, 17002 (2009). 
P. Werner, E. Gull, O. ParcoUet, and A. J. Millis, Phys. 
Rev. B 80, 045120 (2009). 

B. Kyung, S. S. Kancharla, D. Senechal, A.-M. S. Trem- 
blay, M. CiveUi, and G. Kothar, Phys. Rev. B 73, 165114 
(2006). 

D. Senechal and A.-M. S. Tremblay, Phys. Rev. Lett. 92, 
126401 (2004). 

N. S. Vidhyadhiraja, A. Macridin, C. §en, M. Jarrell, and 

M. Ma, Phys. Rev. Lett. 102, 206407 (2009). 

A. Liebsch and N.-H. Tong, Phys. Rev. B 80, 165126 

(2009). 

T. Maier, M. Jarrell, T. Pruschke, and J. Keller, Phys. 
Rev. Lett. 85, 1524 (2000). 
'^'^ S. S. Kancharla, B. Kyung, D. Senechal, M. CiveUi, 
M. Capone, G. Kothar, and A.-M. S. Tremblay, Phys. 
Rev. B 77, 184516 (2008). 

D. Senechal, P.-L. Lavertu, M.-A. Marois, and A.-M. S. 
Tremblay, Phys. Rev. Lett. 94, 156404 (2005). 

''^ M. Capone and G. Kothar, Phys. Rev. B 74, 054513 

(2006). 

M. Aichhorn, E. Arrigoni, M. Potthoff, and W. Hanke, 
Phys. Rev. B 76, 224509 (2007). 
M. CiveUi, Phys. Rev. Lett. 103, 136402 (2009). 
^2 M. CiveUi, Phys. Rev. B 79, 195113 (2009). 

K. Haule and G. Kothar, Phys. Rev. B 76, 104509 (2007). 
O. ParcoUet, G. Biroli, and G. Kothar, Phys. Rev. Lett. 
92, 226402 (2004). 

M. CiveUi, M. Capone, S. S. Kancharla, O. ParcoUet, and 
G. Kotliar, Phys. Rev. Lett. 95, 106402 (2005). 
^® T. D. Stanescu and G. Kothar, Phys. Rev. B 74, 125110 
(2006). 

^'^ S. Sakai, Y. Motome, and M. Imada, Phys. Rev. Lett. 
102, 056404 (2009). 

S. Sakai, Y. Motome, and M. Imada, Phys. Rev. B 82, 
134505 (2010). 

P. Werner, A. Comanac, L. de Medici, M. Troyer, and 
A. J. Millis, Phys. Rev. Lett. 97, 076405 (2006). 
®° P. Werner and A. J. Millis, Phys. Rev. B 74, 155107 

(2006). 

®i K. Haule, Phys. Rev. B 75, 155113 (2007). 

E. GuU, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, 
M. Troyer, and P. Werner (2010), arXiv: 1012.4474. 

63 A. N. Rubtsov (2003), cond-mat/0302228. 

6" E. GuU, O. ParcoUet, P. Werner, and A. J. Millis, Phys. 

Rev. B 80, 245102 (2009). 
®^ M. Ferrero, P. S. Cornaglia, L. D. Leo, O. ParcoUet, 



27 



G. Kotliar, and A. Georges, Europhys. Lett. 85, 57009 

(2009). 

®^ M. Fcrrcro, P. S. Cornaglia, L. Dc Leo, O. ParcoUet, 
G. Kotliar, and A. Georges, Phys. Rev. B 80, 064501 
(2009). 

K. Mikelsons, E. Khatami, D. Galanakis, A. Macridin, 
J. Moreno, and M. Jarrell, Phys. Rev. B 80, 140505(R) 
(2009). 

®* E. Khatami, K. Mikelsons, D. Galanakis, A. Macridin, 

J. Moreno, R. T. Scalettar, and M. Jarrell, Phys. Rev. B 

81, 201101 (2010). 
®^ E. Gull, M. Ferrero, O. ParcoUet, A. Georges, and A. J. 

Millis, Phys. Rev. B 82, 155101 (2010). 
™ G. Sordi, K. Haule, and A.-M. S. Tremblay, Phys. Rev. 

Lett. 104, 226402 (2010). 

A. Macridin, M. Jarrell, and T. Maier, Phys. Rev. B 74, 
085104 (2006). 

J. E. Hirsch and R. M. Eye, Phys. Rev. Lett. 56, 2521 

(1986). 

E. Gull, P. Werner, A. Millis, and M. Troyer, Phys. Rev. 
B 76, 235123 (2007). 

M. Imada, Phys. Rev. B 72, 075113 (2005). 

T. Misawa and M. Imada, Phys. Rev. B 75, 115121 

(2007). 

'■^ M. Sheffer, J. Bascompte, W. A. Brock, V. Brovkin, S. R. 
Carpenter, V. Dakos, H. Held, E. van Nes, M. Rietker, 
and G. Sugihara, Nature 461, 53 (2009). 
J. Joo and V. Oudovenko, Phys. Rev. B 64, 193102 
(2001). 

^® S. Chakraborty, D. Galanakis, and P. PhiUips (2010), 

arXiv:1001.1918. 
'■^ X. Y. Zhang, M. J. Rozcnberg, and G. KotUar, Phys. Rev. 

Lett. 70, 1666 (1993). 
®° P. Werner and A. J. MiUis, Phys. Rev. B 75, 085108 

(2007). 

®i G. Kothar and D. Vollhardt, Physics Today 57, 53 (2004). 

B. Kyung and A.-M. S. Tremblay, Physical Review Letters 
97, 046402 (pages 4) (2006). 

®^ A. H. Ncvidomskyy, C. Schcibcr, D. Scncchal, and A.- 
M. S. Tremblay, Phys. Rev. B 77, 064427 (2008). 
D. B. McWhan, A. Menth, J. P. Remeika, W. F. 
Brinkman, and T. M. Rice, Phys. Rev. B 7, 1920 (1973). 

F. Kagawa, T. Itou, K. Miyagawa, and K. Kanoda, Phys. 
Rev. B 69, 064511 (2004). 

J. H. Shim, K. Haulc, and G. Kothar, Nature 446, 513 
(2007). 

K. Haule and G. Kotliar, Phys. Rev. B 76, 092503 (2007). 
J. M. Ziman, Principles of the Theory of Solids (Cam- 



bridge University Press, Cambridge, 1972). 
R. K. Kaul, A. Kolezhuk, M. Levin, S. Sachdev, and 
T. Senthil, Phys. Rev. B 75, 235122 (2007). 
^° R. Kaul, Y. B. Kim, S. Sachdev, and T. Senthil, Nat Phys 

4, 28 (2008). 

H. Alloul, T. Ohno, and P. Mendels, Phys. Rev. Lett. 63, 
1700 (1989). 

®^ W. W. Warren, R. E. Walstedt, G. F. Brennert, R. J. 
Cava, R. Tycko, R. F. Bell, and G. Dabbagh, Phys. Rev. 
Lett. 62, 1193 (1989). 

G. Aepph, T. E. Mason, S. M. Hayden, H. A. Mook, and 
J. Kulda, Science 278, 1432 (1997). 
■^^ J. L. Tallon, J. W. Loram, G. V. M. WiUiams, J. R. 
Cooper, I. R. Fisher, J. D. Johnson, M. P. Staines, and 

C. Bernhard, Physica Status Solidi (b) 215, 531 (1999). 

5. Sachdev, Quantum Phase Transitions (Cambridge Uni- 
versity Press, New York, 2001). 

A. Kopp and S. Chakravarty, Nat. Phys 1, 53 (2005). 
^'^ H. V. Lohncysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. 

Mod. Phys. 79, 1015 (2007). 
^® S. Roy and A.-M. S. Tremblay, EPL (Europhysics Letters) 

84, 37013 (6pp) (2008). 
®® H. Terletska, J. Vucicevic, D. Tanaskovic, and V. Do- 

brosavljevic (2010), arXiv: 1012.5833. 

M. R. Norman, D. Pines, and C. KaUin, Adv. Phys. 54, 

725 (2005). 

M. R. Norman and C. Pepin, Reports on Progress in 
Physics 66, 1547 (2003). 
"'2 D. Bonn, Nat Phys 2, 159 (2006). 
M. R. Norman, Physics 3, 86 (2010). 

D. Jaksch and P. Zoller, Ann. Phys. 315, 52 (2005). 

M. Lewenstein, A. Sanpera, V. Ahufingcr, B. Damski, 
A. Sen, and U. Sen, Adv. Phys. 56, 243 (2007). 

I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 
80, 885 (2008). 

M. Greiner, O. Mandel, T. EssUnger, T. Hansch, and 
I. Bloch, Nature 415, 39 (2002). 

R. Jordens, N. Strohmaier, K. Gunter, H. Moritz, and 

T. EssUnger, Nature 455, 204 (2008). 

A. Eckardt and M. Lewenstein, Phys. Rev. A 82, 011606 

(2010). 

"° P. Gegenwart, Q. Si, and F. Steghch, Nat Phys 4, 186 

(2008) . 

"1 B. S. Shastry, Reports on Progress in Physics 72, 016501 

(2009) . 

"2 T. Honma and P. H. Hor, Phys. Rev. B 77, 184520 (2008). 



